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ABSTRACT 


Models of noise and dynamic characteristics of a gyro and auto- 
collimator for very small signal levels are presented. Measurements have 
been evaluated using spectral techniques for identifying noise from 
base motion. The techniques have been developed to enable the sensors 
for a relativity experiment to be evaluated in spite of the fact their 
accuracy will exceed any reference standards available. 

Conventional tests of gyroscopes, telescopes, and other precise 
orientation sensors normally rely on more accurate standards to provide 
a reference. These standards are used in an absolute sense, in that all 
errors are attributed to the instruments being tested, including those 
of the calibration equipment. When the precision of the instruments 
under test surpasses that of available laboratory calibration equipment, 
different methods are needed. This situation exists in a test of general 
relativity theory being developed at Stanford University, 

The experiment is constructed to measure the precession, due to 
relativistic effects, of an extremely precise earth-orbiting gyroscope. 

The design goal for non-relativistic gyro drift is 0,001 arcsec per 
year, A telescope, used to maintain the experiment satellite in an 
inertially fixed orientation^ has similar extreme precision. A labora- 
tory test is desired which accounts for environmental effects, especially 
gravity, and isolates inherent sensor errors. 

Developing appropriate test methods using the actual relativity 
sensors was impractical due to the expense and unavailability of proto- 
types. Consequently, an analogous, though less accurate fixed base 
simulator was developed. This simulator was used in developing methods 
of instrument error modeling and performance evaluation applicable to 
the relativity experiment sensors and other precision pointing instruments. 

Two Atlas missile gyros, integrated into a Minuteman 1 guidance 
assembly, and two specially designed autocollimators serve as analogs to 
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the relativity experiment sensors. A three-axis Minuteman I stable plat- 
form provides a movable mount for the sensors. The sensors are used 
to stabilize the platform with less than 1 arcsec jitter. The noise 
equivalent angle of the primary autocollimator is 0.1 arcsec. 

A spectral approach is used in analyzing the fixed base simulator 
signals and developing error models for the sensors, A review of spec- 
tral theory and practical considerations in computing spectral functions 
introduces a discussion of experimental results. Comparisons are made 
of gyro and autocolliraator spectra taken with the instruments both 
resting on quiet test beds and serving as feedback sensors in platform 
stabilization loops. Coherency between signals is used to indicate 
mutual instrument error sources and to provide confidence in the estima- 
tion of system transfer functions. 

Analysis of autocollimator spectra uncovered the presence of a 
platform gimbal resonance. The source of resonance was isolated to 
gimbal bearing elastic restraint properties most apparent at very small 
levels of motion « 100 prad). A model of these properties which 
include both elastic and coulomb friction characteristics is discussed, 
and a describing function developed. 
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I. INTRODUCTION 


A. Problem Statement 

Reliable prelaunch testing and evaluation of extremely precise sat- 
ellite orientation sensors has become an important concern. As the ac- 
curacy of both optical and inertial sensors has improved, it often equals 
or surpasses that of conventional reference instruments. 

In conventional testing, where a precise reference exists, differences 
in output between the test sensor and the reference instruments are usually 
attributed to errors in the test instruments. When a more precise refer- 
ence does not exist, the best source of information about the orientation 
is often the sensor itself. The purpose of this research is to demonstrate 
the effectiveness of a frequency domain approach to orientation sensor 
error modeling and performance evaluation when a precise reference is not 
available. The demonstration is accomplished through experiments on a 
pair of gyros and autocollimators mounted on a stabilized three-axis plat- 
form. 


B. The Stanford Relativity Experiment 

The motivation for this study was a requirement to verify in the lab- 
oratory the design and fabrication of the extremely precise gyros and 
telescope of the Stanford Relativity Experiment. This experiment, first 
proposed in 1960 [SI], is designed to measure the precession of a gyro- 
scope predicted by the general theory of relativity. The accuracy will 
be sufficient to distinguish between Einstein's theory and other theories 
which have been proposed. 

Precession is measured relative to a telescope pointed at a star. 

As outlined in [Bll], two general relativistic precession effects are of 
interest. The first, known as geodetic precession, is due to motion of 
the gyros through a gravitational field. For a 500 nautical mile altitude 
earth polar orbit, the predicted precession is 6.9 arcsec per year. An 
even smaller effect, called motional drift, is due to rotation of the 
gravity field. The secular portion of this effect for the same orbit is 
0.05 arcsec per year. This experiment differs from all other proposed 
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relativity experiments in that it is the only one which can measure the 
effects of a gravity field produced by a rotating body. In order to ob- 
serve the latter precession to 2$, accuracy, a design goal of 0.001 arcsec 
per year of drift due to Newtonian torques has been established. 

In addition to the extreme accuracy required of the gyros, the null 
stability of the telescope must also be held to 0.001 arcsec. To meet 
this requirement, the telescope is limited to a linear dynamic range of 
0.1 arcsec and is of special construction. 

Achieving such accuracies, more than six orders of magnitude better 
than existing instruments in the case of the gyros, obviously requires 
special techniques. The two principal aids are the low g orbital envi- 
ronment and use of liquid helium temperatures. The environment is pro- 
vided by placing the sensor package inside a dewar filled with superfluid 
helium. At these temperatures creep is virtually nonexistent, providing 
excellent mechanical stability, and the superfluid helium is a perfect 
heat conductor which minimizes temperature gradients. A superconducting 
niobium shield around the experiment package eliminates magnetic inter- 
action torques between the gyro and the earth's magnetic field. Super- 
conducting electronics are also used in the essentially torque-free gyro 
pickoff, and helium boiloff is used in linear differential thrusters [B12] 
to provide attitude control to within 0.05 arcsec rms. 

The gyro rotors are made of fused quartz covered with thin metallic 
superconducting layers. The rotors are electrostatically suspended with 

g 

a nominal orbital preload of 1 x 10 g to prevent satellite disturbances 
from causing contact with the gyro. Gravity-gradient drift torques are 
reduced to acceptable levels by producing the rotors spherical and homo- 
geneous to within 10 . Details of the design and mission can be found 

in [E2], [Bl], and [Bll]. 

C. Fixed Base Simulator 

Although the precision just indicated for an orbital environment 
cannot be obtained in the laboratory, primarily due to the one-g environ- 
ment, it is nevertheless important to have measures of performance to 
ensure that the instruments have been properly fabricated and assembled. 
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Proper operation of the instruments requires that they be in a liquid 
helium dewar, where they have poor accessibility. It may not be possible 
to keep the instruments within their linear range due to mechanical shifts 
and thermal distortion without some means of sensor output feedback. 
Further, since drift due to gravity will be the dominant effect seen in 
the laboratory, it is desirable to be able to change the orientation of 
the sensor package relative to the local vertical in order to separate 
various error terms. 

For these reasons, a fixed base simulation using a gimballed dewar 
mount and an artificial star assembly to provide a reference target for 
the telescope is being constructed to test the experiment sensors. 

The complexity, inflexibility, and expense of working with the Rela- 
tivity Experiment package in its helium environment made a preliminary 
study using similar but less sensitive hardware desirable. Consequently, 
a relatively low cost and considerably more flexible analogous fixed base 
simulator was constructed. It contains a pair of Atlas missile two-degree 
of-freedom floated gyros which provide an analog to the Relativity Experi- 
ment free-rotor gyros. A small specially built autocollimator and refer- 
ence mirror are similar in behavior to the telescope and artificial star 
of the main experiment. The dewar mount is replaced in the fixed base 
simulator by a three-axis gimballed platform on which are mounted the 
gyros and autocollimator. A detailed description of the simulator hard- 
ware is given in Chapter II. It has been used to develop and demonstrate 
test methods which can be applied to the Relativity Experiment and other 
fine-pointing experiments. 

D. Related Work 

1. Isolation Platform 

The problem of providing a quiet, stable test bed to evaluate 
precise motion sensors has received considerable attention. Concrete 
piers in seismically quiet areas, passive isolation using pneumatic or 
mechanical spring supports, and active servo-controlled test platforms 
have all been tried. 

The Martin Company [Ml] has constructed a large (23 ft x 39 ft 
X 8 ft) concrete block resting on bedrock for testing inertial guidance 
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components. Its location, just outside Denver, Colorado, was selected 
after extensive seismic studies because of its relative quietness. Flack 
and Dickie [FI] studied the effects of ground noise on accelerometer and 
gyro testing when the Singer- Kearfott Company decided to construct a 
dynamic test facility in close proximity to an existing inertial test 
facility. Of particular interest was the noise transmission characteris- 
tics between the facilities. 

Numerous passively sprung isolation platforms have also been 
constructed. These platforms are effective in attenuating high frequency 
disturbances, but move with the underlying structure at low frequencies. 

A six-degree-of-freedom isolation platform which is suspended by a pen- 
dulum from an isolating suspension system has been used for precision 
testing at the Aerospace Guidance and Metrology Center at Heath, Ohio 
[P3] . The Seiler Laboratory at the Air Force Academy has also developed 
an isolation platform, described by Wittry [W6], which is supported by 
twenty pneumatic isolators. Tsutsumi [T2] reviews these and other isola- 
tion platforms in a geokinetics review paper. 

When greater precision and low frequency isolation has been 
required, actively servo-controlled platforms have been developed. 

Tsutsumi [T3] developed a ground tilt isolation platform which used 
electronic tilt sensors and hydraulic actuators to provide isolation 
from low frequency floor motion. Weinstock [Wl] improved this system 
by using a gyro for higher frequency sensing and control with a tilt sen- 
sor eliminating gyro drift errors. DeBra, Mathieson and Van Patten [D2] 
developed a two-degree-of-f reedom table leveling system employing direct 
electromagnetic force actuators for low frequency attenutation. These aug 
ment the primary support and high frequency attenuation afforded by 
Belleville type springs. A temperature controlled level senses the local 
vertical. Lorenzini [L2] describes improvements to the Seiler Laboratory 
platform which control the platform actively in all six degrees-of-f ree- 
dom. A combination of pneumatic cylinders and electromagnetic shakers 
provide control forces, while tiltmeters, angular accelerometers, seismo- 
meters, and a gyrocompass sense the platform motion. Another tilt isola- 
tion platform has been constructed at the Newark Air Force Station, Ohio 
[P2]. 
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2 . 


Simulation 




The primary purpose of all isolation platforms is to reduce 
instrument output signals caused by test base motion to levels substan- 
tially lower than those created by instrument errors. Other applications 
require accurately known base motion in order to calibrate the motion 
sensing instruments and to perform system checkouts. Such fixed-base and 
moving-base simulations have become an integral part of the development 
of nearly all aerospace vehicle and transportation systems. Reinel [Rl] 
describes a typical satellite fixed base simulation used in testing the 
European Symphonic satellite. The sensors were mounted on a gimballed 
platform whose orientation is continuously controlled by servo loops 
in response to angle and angular velocity or rate commands from an analog 
computer. Platform angles are sensed by potentiometers or encoders while 
rates are sensed by tachometers. Differences between these signals and 
the commands are nulled by accelerating the platform gimbals with electro- 
hydraulic actuators and an electric motor. 

Davis [C8] describes a precision satellite attitude determina- 
tion system called SPARS which uses strapdown gyros for short term con- 
trol and discrete star sightings to correct for gyro drift. The laboratory 
test of this combined optical-inertial system is described by Hvoschinsky 
and Horiuchi [H2]. A sophisticated six-degree-of-f reedom hybrid simula- 
tion is also discussed by Hogan and Welch [Hi]. A salient feature of all 
these simulations is that primary orientation sensing is provided by test 
fixture transducers and not by the instruments under test. This feature 
is not present in the Fixed Base Simulator described in this dissertation. 

In very accurate pointing systems, the components under test 
may far exceed the accuracy of available test fixture transducers. Hence 
there is no absolute laboratory standard for comparison and calibration. 
Since typical operating signals of these systems are less than 1 arcsec, 
it is necessary to do modeling and testing with signals of comparable 
amplitude to obtain realistic results. Even though larger signals are 
easier to measure accurately, the increased motion can alter significantly 
the physical characteristics of the hardware components. 

A fairly recent branch of the simulation field known as "model- 
following" is concerned with accurate simulation of systems by means of 
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other possibly dissimilar systems. Tyler [T4] discusses an optimal con- 
trol approach to the problem. Other interesting aspects of this area 
have been studied by Erzberger [El], Asseo [Al], Curran [C6], and Winsor 
and Roy [W5] . 

3. Data Analysis 

Analysis of orientation sensor signals using spectral techniques 
to determine instrument errors is quite recent although the theory has been 
available for several decades [W4], [Jl], [Nl]. Earlier efforts to do so 

were limited by the expense of computing the spectra, a situation which 
has been remedied by the development of high speed computers and efficient 
algorithms. An interesting application of the modem digital spectral 
techniques used in this dissertation, is the analysis of the voices of 
killer whales discussed by Singleton and Poulter [S2]. Their paper gives 
an extensive discussion of the methods used in obtaining the voice spectra. 
Continuous recordings taken from underwater microphones were broken into 
short segments and the spectra computed. By comparing the spectra, speech 
patterns could be discerned. To name a few other applications, spectral 
techniques have been used in the study of electrocardiograms, ocean waves, 
atmospheric radio propagation, and seismic events. 

More closely related to the present study is an analysis of the 
effects of vibration on gyro testing performed by Crowly and Ossing [C5] . 
Vibration of the test platform was measured by pendulous seismometers 
and the sensed angles subtracted from the gyro output. The coherence 
between the gyro and sesimometer signals indicated the relative contribu- 
tion of vibration to the test gyro signal. 

Another closely related experimental study of the spectral char- 
acteristics of gyro drift was conducted by Truncale, Koenigsberg and Harris 
[Tl]. A gyro with low mid-frequency noise (0.01 to 100 Hz) was desired 
for use in the Orbiting Astronomical Observatory. Gyros from several dif- 
ferent manufacturers were tested. Analyzing the frequency distribution 
of drift noise proved powerful in selecting a suitable gyro. A source of 
concern in the tests which was discussed extensively was the uncertainty 
in the amount of output due to base motion. 
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Although many other spectral studies of inertial instruments 
are known to have been performed, there remains a paucity of experimental 
results available in the open literature particularly for small signal 
levels. 


E. Summary of Research 


The principal contributions of this research are 

• the construction of a fixed base simulator to aid in developing 
methods of testing precise optical and inertial orientation 
sensors. 

• fabrication and design of a precise autocollimator 
linearity better than 0.1$ and noise equivalent angle of 
0.5 [irad. 

• the analysis of spectra of the precise autocollimator and an 
Atlas missile two-degree-of-f reedom gyro and the development 
of noise models for each. 

• the measurement of small motion behavior of ball bearings and 
the development of a describing function. 

. the demonstration of an approach to develop error models of 
optical and inertial sensors. 

• the establishment of a method for evaluating the Relativity 
Experiment instruments. Programs for numerically carrying 
this out have been written and verified using this fixed base 
simulator. 
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Chapter II 


FIXED BASE SIMULATOR 

In this chapter, the Fixed Base Simulator developed as an analog 
to the Stanford Relativity Experiment laboratory test is described. An 
overview of the system is first presented. Since the configuration of 
the simulator was heavily influenced by hardware constraints, a discus- 
sion of its development follows. The subsystems introduced in the de- 
velopment section are then discussed in detail. Finally, the modes in 
which the simulator can be operated are described. 


A. Overview 

The Fixed Base Simulator consists of a three-axis stable platform 
assembly, motion sensing instruments, and electronics, Fig. 1. The plat- 
form is a modified government surplus Minuteman I stable platform on which 
are mounted two Atlas missile gyros and a specially designed two axis auto 
collimator. The platform and an autocollimator reference mirror are in- 
stalled on a 24 in. rotary table, which serves as an earth orientation 
reference. A second miniature autocollimator, necessary to give the third 
optical measurement of orientation, is attached to the rotary table and 
uses a reference mirror mounted on the gimballed platform. The platform 
and sensors are shown in Fig. 2. 

Minuteman I electronics used in the simulator are located in their 
original guidance assembly chassis. Fig. 3, while newly developed elec- 
tronics are integrated into a separate control console. Since the gyros 
and autocollimators are not Minuteman I equipment, extensive electronic 
and mechanical modification and new development was required to make them 
compatible with the gimballed platform. Original Minuteman electronics 
used are (1) various power supplies, (2) gyro pickoff preamplifiers, (3) 
pickoff excitation amplifier, (4) pitch transformation resolver, and (5) 
control amplifiers, consisting of demodulators, compensation networks , and 
power amplifiers. New equipment consists of (1) gyro torquer and caging 
amplifiers, (2) heaters, (3) preamplifiers, and (4) torquer magnet excita 
tion; (5) autocollimator circuitry, (6) several dc and ac precision power 
supplies, (7) mode switching logic and relays, (8) interconnect cabling, 
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FIG. 1 FIXED BASE SIMULATOR. The control console and Minuteman I 
electronics support and control the platform and sensors. 
Analog recorded data is digitized before processing. 
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FIG. 2 PLATFORM AND SENSORS ASSEMBLY. The two autocollimators 
or the two gyros provide the reference required for 
stabilization. 
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FIG. 3 MINUTEMAN I ELECTRONICS TOROID. Several modules are 
used to provide power and control platform motion. 
The D-17 computer is not used. 


and (9) instrumentation. These are all discussed in greater detail later 
in this chapter. 

In all modes of operation, the platform is stabilized using either 
the gyro or the autocollimator signals. Figure 4 shows the basic oper- 
ation of the Fixed Base Simulator. 



Fig. 4. FIXED BASE SIMULATOR FUNCTIONAL DIAGRAM. 


Motion of the platform is sensed by the gyros and the autocollima- 
tors, producing ac signals proportional to angular displacement. The 
instrument signals are amplified and demodulated and then sent to the 
instrumentation equipment. Also, depending on the mode of operation, 
either the gyro signals or the autocollimator signals are passed through 
compensation networks and amplified to drive the platform gimbal torque 
motors. The direction of torque is that required to null the controlling 
sensor output. Simultaneously, signals from the noncontrolling sensors 
are used as inputs to the gyro caging amplifiers, required to keep the 
gyro pickoff in the small linear region about null, or to maintain opti- 
cal reference. 
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Command inputs are introduced in three different ways. First, the 
outputs of the autocollimators can be biased directly. Second, bias 
torque levels can be applied to the gyros through the torquer amplifiers. 
Finally, command torque signals can be added to the compensated error 
signals at the inputs to the gimbal torquer power amplifiers. 

Several modes of operation allow selection of the controlling sen- 
sors, on and off switching of gyro spin motor power, and inhibiting of 
caging amplifier outputs. 


B. Simulator Development 

Before giving a more detailed discussion of the system, the evolu- 
tion and development of the Fixed Base Simulator is described both as 
a matter of interest and to give a better understanding of the final 
configuration. This description is basically chronological in order, 
and attempts to outline the rationale and constraints which affected the 
course of development. 

1. Platform Selection 

The first basic requirement of a fixed base simulator was a 
movable platform to which motion and attitude sensors could be mounted. 
To speed development, existing hardware was sought. It was found that 
several aircraft and missile stable platforms were available either as 
government surplus or commercially. Of these platforms, the ten -year- 
old Minuteman I was selected for several reasons. The primary advantage 
of this system is its "inside-out” gimbal arrangement, in which the mo- 
tion sensors are physically located outside of the gimbal bearings and 
no gimbal rings are required. This arrangement allows mounting of test 
instruments to the control element without volume constraints. The di- 
rect drive dc gimbal motors do not inherently limit pointing accuracy, 
and the precision gimbal bearings do not present a serious stiction 
problem. Furthermore, two complete guidance electronics packages and 
several stable platforms were available as surplus equipment, eliminat- 
ing any spare parts concern. 
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The major disadvantage of the Minuteman I system was the un- 


availability of the gyros and velocity meters due to their security clas- 
sification. 

The Minuteman I two-degree- of -freedom air bearing gyros would 
have been ideally suited for use with the simulator, but it was felt un- 
likely that the instruments would be declassified in time to be of use. 

It was further felt that other gyros could be made compatible 
with the Minuteman system, and consequently, an order was placed for the 
equipment . 


2. Sensor Selection 

The second basic requirement of the simulator was sensors 
closely analogous to those of the Relativity Experiment; i.e., free rotor 
gyros and a telescope. Consideration was first given to the gyros. 

a. Gyro 

A search for gyros revealed no suitably priced free rotor, 
or gimballed two-degree-of- freedom gyros other than several fifteen year 
vintage Atlas missile instrument gyros already possessed by the Stanford 
Guidance and Control Laboratory. The gyros are large (4.6 kg), two-de- 
gree-of -freedom, floated instruments. 

Available documentation on these instruments consisted 
solely of a brief familiarization manual. No wiring diagrams, mechanical 
drawings, detailed performance specifications, or operating manual were 
found. Nevertheless, since the gyros were immediately available, and 
were reasonably analogous to the Relativity Experiment gyros, it was de- 
cided to attempt to integrate them into the Minuteman I system. 


b. Autocollimator 

Before proceeding with the gyro modification, however, the 
question of a suitable analog to the Relativity Experiment telescope had 
to be resolved. Use of a telescope would have required development of an 
artificial star, and made the system subject to angular errors arising 
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from translations of the telescope and artificial star. An autocollima- 
tor, which provides its own light reference, has neither of these disad- 
vantages, yet is sufficiently similar to a telescope to make meaningful 
studies. A survey showed that no autocollimator meeting the simulator 
requirements was commercially available. Since an in-house capability 
existed to design and build a suitable model, it was decided to develop 
one rather than compromise simulator design goals. Design and fabrica- 
tion of the instrument was postponed until the rest of the system was 
near completion, since no major development problems were anticipated. 


3. Minuteman I Guidance Assembly 

Delivery of the Minuteman I guidance assemblies ushered in an 
extremely interesting and enlightening phase of the system development. 
The assemblies arrived just as they were taken from the missile, with 
the exception that the inertial instrument, fan, battery, and ordnance 
switching unit had been removed and the hermetic seal broken. 

Each cylindrical assembly is covered by a reflective surface 
to prevent excessive thermal radiation to internal parts and to provide 
a hermetic seal. Inside the cover is a hollow, twelve-sided, cylindrical 
structure, referred to as the toroid, in which the guidance computer and 
stable platform electronics are located. The stable platform rests on 
a heat exchanger surrounded by the toroid. The toroid is separable into 
two parts: the D17 guidance computer, and the platform electronics. Al- 

though functional, the computer is not employed in the simulator because 
of the amount of supporting equipment and software which would have been 
required . 

Mounted on the stable platform is an alignment unit consisting 
of bubble levels for vertical reference and a mirror for use as an opti- 
cal azimuth reference. 

A period of hardware familiarization followed, in which the 
physical layout, external connectors, and internal wiring were studied. 

Before power could be applied to the system, cooling require- 
ments had to be evaluated. The original Minuteman I ground cooling sys- 
tem was complex and required pumping refrigerated water through the 
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assembly in addition to circulating pressurized helium by means of a 
vaneaxial fan. 

Laboratory cooling requirements were judged not as severe 
since a hermetic seal is not necessary. However, the amount of cooling 
required was difficult to analyze. Furthermore, a decision on whether 
or not to use the guidance computer had not yet been made. To ensure 
adequate cooling, design of an air blower ducting system was therefore 
conducted in parallel with other tasks. 

The toroid was removed from the missile mounting structure and 
bolted to a circular plywood frame equipped with casters to provide mo- 
bility. A blower was attached to the "bonnet” on top of the toroid which 
originally held the vaneaxial fan. Sheet metal skirts were placed around 
the toroid base to cover the gap between the bottom of the toroid and the 
plywood frame. In this configuration, air was forced down through the 
center of the hollow toroid, deflected by the plywood and sheet metal 
skirts, and circulated back up past the electronics modules mounted in 
the vented toroid walls. Later, after deciding not to use the computer, 
it was found that operating temperatures were acceptable without using 
the cooling system. As a consequence, the blower was removed from the 
toroid in the final configuration. 


4 . Power 

Concurrent with the cooling system development, a power source 
for the platform electronics was sought. The Minuteman system requires 
28 VDC, but current requirements without the computer were not accurately 
known. Estimates ranged from 10 to 30 A. A reasonably priced ($25) 
power supply meeting the maximum requirements was procured. The unit is 
a well-regulated, safeguarded, Autonetics surplus supply which had orig- 
inally been a part of the Minuteman I ground equipment. The sole draw- 
back of this power supply was that it required 115 V, 30 , 400 Hz input 
power which was. not available. To provide 400 Hz power, a pair of 60 Hz 
to 400 Hz motor- generator sets, U. S. Marine Corps surplus, were reclaimed 
from Stanford University storage. The required input power to the motor- 
generators, 115 V, 30 , 60 Hz, rather than the 208 V, 30 power available 
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in the laboratory, was provided by three large power transformers again 
found in storage. Connecting these transformers in a Wye configuration 
provided the necessary voltage stepdown. All of the components were 
cabled together and satisfactory power was established. 

With adequate power available, studies were made of the plat- 
form electronics voltage levels, waveforms, and other characteristics. 
These studies were necessary since only module schematics and no detailed 
performance specifications were available. The primary objective was to 
determine the interface with the rest of the simulator. 

5. Gyro Operation 

At the same time, work began on determining the operating char- 
acteristics of the Atlas gyros. The work was severely hampered by the 
lack of documentation. It is imperative to note that the professional 
assistance of Dr. W. G. Pondrom was crucial in the successful completion 
of this task. 

After extracting two gyros from an Atlas stable platform, it 
was first necessary to determine the function of each connecting wire. 

A gyro wiring schematic was eventually deduced from ohmic measurements. 

A gyro test breadboard was then constructed which allowed activation of 
the pickoffs and torquers and application of power to the spin motor. 


a. Pickoff and Torquer Coupling 

It was soon deduced that the pickoff and torquer coils 
were imbedded in a common magnetic structure. This implies that non-can- 
celling ac currents in the pickoff coils, which occur when the pickoff 
is off-null, will induce currents in the torquer coil, thereby loading 
the pickoff circuit. The result is poor off-null pickoff performance. 

In the intended application, this characteristic was not important since 
the gyro was held at null by the servo control loop. 

However, simulator requirements demanded good off-null 
performance, so equivalent current source torquers were needed to elimi- 
nate pickoff loading. Design of torquer amplifiers was postponed while 
other characteristics were determined. 
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b. 


Pickoff Modification 


It was necessary to make the gyro pickoffs compatible 
with the Minuteraan electronics. Atlas pickoff excitation frequency is 
400 Hz while the Minuteman I frequency is 4800 Hz. In addition, the 
Atlas excitation requires 4.65 VDC power for the torquer electromagnet 
with one side grounded, while the Minuteman gyro has no dc excitation 
and is capacitively coupled. Finally, Minuteman excitation amplifier 
voltage is 100 V rms , while the Atlas originally used 8 V. (The 8 V 
amplitude was later changed to 30 V to increase pickoff gain.) Exten- 
sive modifications were clearly required. Detailed documentation of 
these modifications is presented in the subsystems section later in this 
chapter. 


c. Spin Motor Power 

The next design step was to provide spin motor power to 
the gyro. Available Minuteman power was 400 Hz, 22 V, 30 quasi-square 
wave (+,0,-,0, +,...). Required Atlas power is 24 V, 30 sine wave with 
one side grounded. Isolation and voltage stepup to achieve the appro- 
priate rms power were accomplished by hand winding three bifilar coils 
on ferrite cores. These transformers were placed in a custom box and 
potted in beeswax. 


d . Breadboard Test 

At this point the gyro test breadboard was sufficiently 
complete to fully operate the gyro, with the exception of the heater. 

The Minuteman caging amplifier, used in controlling the redundant fourth 
gyro axis, was taken from each of the two guidance assemblies. The gyro 
pickoffs were connected to the amplifier inputs and the torquers to the 
outputs through isolation transformers. The caging amplifier outputs 
are very nearly current sources so pickoff loading was not a problem. 

The gyros were then successfully operated closed loop, with the pickoff 
outputs held at null. Keeping one gyro axis caged, rough estimates of 
the pickoff and torquer scale factors were made using earth rate as a 
reference. These estimates were used in the torquer amplifier design 
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and also indicated the need for pickoff preamplifiers to achieve suffi- 
cient platform stabilization loop gain. 

Enough information had been gained at this time to proceed 
with the electronics design. Mechanical design of two adapter brackets 
to hold the gyros on the platform was also begun. 


6. System Configuration 

While this work was going on, system layout design was initi- 
ated. Connector pin assignment, cabling requirements, grounding, shield- 
ing, switching, and monitoring were each planned in detail prior to sys- 
tem assembly. 


7. Autocollimator Design 

As system layout and electronics design proceeded, and after 
the gyro mounts had been completed, design of an autocollimator was 
started. Once the optics and light source had been selected, detailed 
mechanical drawings were made. An adjustable reference mirror mount 
also was designed. 


8. Rotary Table Platform Mount 

Questions about where to mount the platform and how to orient 
it to locally vertical and polar positions were resolved when a 24 in. 
rotary table with base was made available by the Stanford Physics Depart- 
ment. A steel plate with bolt brackets was fabricated to adjust the 
table base in azimuth and to level the table. The table and plate were 
positioned at a point where theodolite sightings of Polaris could be made 
through a special window in the laboratory ceiling. 


9. Platform Stabilization 

Upon completion of the electronics design and system wiring, 
but before the autocollimator was completed, the platform was stabilized 
using the gyros. At first, the system was only marginally stable due to 
the decrease in loop gain caused by the increased moments of inertia of 
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the modified platform. Stability of the system was improved by increas- 
ing the gain of the pickoff preamplifiers from 10 to 40. Bandwidth of 
the servo loops was observed to be approximately 15 Hz. 

The autocollimator was completed several weeks after the ini- 
tial gyro stabilization. This autocollimator and a second unit designed 
on another Stanford project were incorporated into the system. Initial 
attempts to stabilize the platform using autocollimator signals failed 
due to instability induced by flexibility in the primary autocollimator 
mount. The instability was eliminated by filtering, reducing the band- 
width of the autocollimator signal below the dominant bending frequency. 

10. System Verification 

The temporary autocollimator reference mirror mount which had 
been used up to this point was replaced by a newly designed unit. The 
new assembly was attached to the rotary table and aligned. After align- 
ing the mirror and bringing autocollimator and gyro axes into coincidence, 
basic data on scale factors, drift rates, and noise levels were taken. 

A peculiar anomaly appeared in the gyro and auto collimator out- 
put signals upon close inspection of the data. It was discovered that 
large, very low frequency disturbances were appearing in the instrument 
signals. In particular, when the platform was stabilized using either 
the autocollimator or gyro signals, the offset in the noncontrolling in- 
strument varied approximately as a square wave with 20 arcsec amplitudes 
and 40 to 400 sec periods. 

The cause of the disturbance was eventually isolated to an in- 
teraction between the 400 Hz quasi-square wave power source for gyro spin 
and the 1.2 MHz crystal oscillator, used in modulating power to the auto- 
collimator light source. The two separate ac sources would lock together 
for many seconds and then slip relative to each other. When this happened 
there was a sudden large change in offset in the autocollimator output. 

A variable frequency generator was substituted for the 1.2 MHz 
crystal and was varied from 0.9 to 1.5 MHz. The noise output of the in- 
strument increased rapidly as the drive frequency approached 1.2 MHz in 
either direction, but just as it reached 1.2 MHz, the random noise disap- 
peared and the low frequency square wave offsets appeared. The 1.2 MHz 
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crystal oscillator, which had been selected at random, had the ideally 
wrong operating frequency. By varying the frequency, a noise minimum 
was found at 1.32 MHz, and a new crystal at this frequency was installed. 

New measurements were made of gyro bias, torquer scale factor, 
and gyro pickoff scale factor. The bias and torquer scale factor were 
deduced by turning the rotary table about a vertical axis while under 
autocollirnator control, and measuring the torquer amplifier voltages 
required to keep the gyro pickoff s at null. Gyro pickoff scale factors 
were found by putting known step changes into the torquers and measuring 
the time required for the pickoff voltage to change a known amount. Ex- 
cellent consistency of results indicated the system was operating prop- 
erly. 


C. Subsystems 

All of the major components of the Fixed Base Simulator have been 
introduced above in the discussion of system development. Since an in- 
timate knowledge of most of these components is necessary to interpret 
the results presented later in this dissertation, details on their con- 
struction and operation are now provided. Circuit schematics of newly 
developed electronic modules are included in the Appendices. 

1 . Sensors 

The sensor subsystem consists of the two gyros, and the primary 
and azimuth orientation autocollimators. 


a. Gyros 

The ARMA Atlas missile gyros used in the simulator, and 

pictured in Fig. 5, are two-degree-of-f reedom, floated inertial grade in- 
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struments. Angular momentum is 10 X 10 gm-cm /sec. 

A spherical float is centered inside the housing by means 
of two orthogonal sets of flexible wires, which also serve as power leads. 
One set connects the float to an intermediate gimbal ring; the other set, 
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FIG. 5 ARMA ATLAS MISSILE TWO-DEGREE -OF-FREEDOM GYROSCOPE. Viscosity 
of the flotation fluid creates low frequency cross-axis 
coupling. Except for this, the gyroscopes are closely analogous 
to the free rotor gyros of the Stanford Relativity experiment. 



oriented 90° away, connects the gimbal ring to the housing. A rotational 
degree of freedom exists about each set of wires. Within the float, a 
pair of flywheels are driven by two 400 Hz, 4 pole induction motors. Also 
contained in the float are a pair of electromagnets, one at each end of 
the spin axis, which are excited by both ac and dc voltages. A pair of 
coils are mounted on the housing directly across from the electromagnets 
to provide two axis pickoff and torquing functions. Direct current 
through the torquer coils sets up a magnetic field which reacts with the 
dc electromagnet field to produce precessional torques. AC excitation 
of the electromagnets induces currents in the pickoff coils, wound in 
the form of figure eights. When the float is off null, these coils de- 
velop voltages proportional to the relative angular displacement of the 
float. 

Heating coils and temperature sensing resistors are mounted 
on the housing to keep the flotation fluid at the neutrally buoyant tem- 
perature. The Atlas platform was temperature controlled as a unit so 
there were no gyro heater electronics. Thus, new temperature control 
circuits were designed. 

Measurements were made of the viscous friction and elastic 
restraint parameters. This was accomplished by torquing the float with 
the wheel turned off. Observing changes in steady state pickoff voltage 
with known small changes in torque yielded estimates of the elastic re- 
straint coefficient h/k. Measured values are 424, 708, 1072, and 515 
sec for the X, Y, Z, and X (redundant) gyro axes, respectively. 

KD 

Similarly, observing the time rate of change of pickoff voltage with large 
changes in torque, gave measurements for the dimensionless viscosity co- 
efficient h/c. These are 1.96, 21.0, 35.7, and 2.43. 

(1) Torquer Amplifiers . As discussed in the previous 
section, new torquer electronics had to be designed to overcome undesir- 
able coupling between the pickoff and torquer. A broadband current 
source was initially contemplated to eliminate the coupling. However, 
after experiencing difficulty with amplifier oscillation, the amplifier 
was designed to have very high output impedance only at the critical 
4800 Hz pickoff excitation frequency and moderate impedance elsewhere. 
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Each torquer amplifier consists of a high bandwidth, low drift, chopper- 
stabilized preamplifier feeding a high-current driver with a twin-T cir- 
cuit in the driver's feedback path. The twin-T provides high loop gain 
at 4800 Hz. A current sensing resistor in series with the load provides 
the necessary feedback to make the combined circuit a high bandwidth (70 
KHz) current amplifier with excellent performance at 4800 Hz. 

The torquers provide extremely stable torquing current 
with drift rates less than 1 (iA/day, thereby eliminating the torquer am- 
plifiers as a significant source of gyro drift. They also have a 140, 000: 1 
linear dynamic range with a maximum torquing capability of 3.5 times earth 
rate. This large dynamic range allows compensation of large biases while 
not masking the more interesting small changes in required torquing levels. 
Torquer amplifier output can be precisely biased by adjusting a ten turn 
potentiometer connected to a precision voltage regulator. The wiper volt- 
age provides input to the torquer amplifier. Closing a control panel 
switch allows maximum torquer output for gyro slewing. Slewing capabil- 
ity is required to quickly reorient the platform without forcing the gy- 
ros into their stops. 

When all four torquer amplifiers were initially connected 
to the same power supply, heavy oscillation resulted. The oscillation, 
caused by feedback through the small power supply impedance, was elimi- 
nated using four decoupling networks. A low pass input filter was also 
required to keep the Minuteman redundant axis caging amplifier from sat- 
urating. 

Each amplifier has two inputs: one driven by a caging 

amplifier, the other by the biasing circuit. The two scale factors , 

v 

K T g are -6.04 and -10.00 mA/V, respectively. Another scale factor, 
which has proven to be more convenient is use, relates the induced gyro 
precession rate to the voltage across the 92.1 ohm torquer amplifier 
current sensing resistor. This parameter was deduced from a series of 
measurements of the torque voltage, taken at 30° azimuth increments with 
the spin axis vertical, required to hold the pickoff at null. Fixed and 
gravity sensitive drift effects are held constant in this way, leaving 
only the known horizontal component of earth rate in the gyro input axes 
as a variable. The values of this parameter, K_ , are 31.6, 33.4, 34.1, 
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32.0 (irad/sec-V for the X, Y, Z, and gyro channels, respectively. 

All other desired torquer scale factors are easily derived using the 
above values. 


(2) Pickoff Excitation and Preamplifiers . In the discus- 
sion on system development earlier in this chapter, several problems in 
making the Atlas gyro pickoffs compatible with the Minuteman electronics 
were pointed out. The Atlas gyro requirements for both ac and dc electro- 
magnet excitation, the phase shift due to a higher excitation frequency, 
and a lower excitation level were the primary concerns. New requirements 
of the pickoff, not related to incompatibility, were to reject noise and 
to establish the proper loop gain for system stability. Electronic modi- 
fications separated into two areas: alteration of the excitation signal 

and the development of a preamplifier. Each is discussed in turn. 

The 100 V rms Minuteman pickoff excitation amplifier out- 
put is used to drive the primaries of two stepdown isolation transformers. 
One end of each secondary is connected to a 5 VDC power supply and the 
other end to a gyro excitation coil. This configuration simultaneously 
isolates the torquer magnet dc excitation from the amplifier and reduces 
the excitation voltage to a nonsaturating level. The amplifier is mounted 
at the base of the platform to isolate its high voltages from other elec- 
tronics. The transformers are conveniently mounted on one of the gyro 
adaptors on the stable platform. 

Compensation of the 40° phase lag could have been accom- 
plished either in the excitation or in the pickoff output. However, un- 
certainties in the magnitude and stability of the Q of the excitation 
isolation transformer made it desirable to shift the phase at the input 
to the preamplifier. Each preamplifier, therefore, consists of a narrow 
passband filter centered at the 4800 Hz excitation frequency, followed by 
a unity gain phase lead network at the input to an amplifier with single 
ended input and differential output. 

The gain of the three primary preamplifiers is 40, while 
the redundant channel is 10. An isolation transformer was also required 
at the output of the redundant channel preamplifier to make it compatible 
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with the Minuteman caging amplifier. To minimize noise pickup, the pre- 
amplifiers are located in a box atop one of the gyro brackets. 

(3) Pickoff Postamplifiers . The amplified signals from 
the new pickoff preamplifiers are sent to the Minuteman electronics tor- 
oid via the control console. Inside the console the three primary chan- 
nel signals are tapped so that each signal is sent to two separate Min- 
uteman ac amplifiers. These amplifiers, which formerly were the gyro 
and velocity meter (VM) pickoff amplifiers, were renamed postamplifiers 
to distinguish them from the newly developed preamplifiers. The VM and 
gyro postamplifiers are all identical ac differential amplifiers with a 
gain of 220 and negligible phase shift. 


(4) Pitch Transformation Resolver . Signals sent to the 
VM postamplif iers are amplified and immediately demodulated in the Min- 
uteman VM control amplifiers. However, two of the three signals sent to 
the gyro postamplifiers pass through the platform pitch transformation 
resolver before being demodulated. The resolver couples the signals by 
the relationships 


x' = x cos 9 + y sin 9 
y ' = -x sin 9 + y cos 9 

where 9 is the inner platform gimbal angle. When the angle 9 is 
zero, both the VM and gyro demodulated signals are the same except for 
electronic errors. The transformed signals are required for platform 
stabilization when the gimbal axes are not orthogonal. 


(5) Pickoff Demodulators . The outputs of the postampli- 
fiers and the pitch transformation resolver are finally sent to the con- 
trol amplifiers for demodulation. A 4800 Hz square wave signal, synchro- 
nous with the signal used by the pickoff excitation amplifier, provides 
the reference. Formerly, the demodulated signals went directly to the 
platform and VM servo compensation networks. In the Fixed Base Simulator, 
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however, the signals are sent to instrumentation jacks and to relays. 
The relays allow the signals to be opened selectively, transferred to 
the caging amplifiers, or sent back to the control amplifier compensa- 
tion networks. 


(6) Pickoff Scale Factors . Measurements of the overall 
pickoff scale factors are based on knowledge of the torquer scale factor, 
which in turn is based on earth rate as the primary reference. The pro- 
cedure is as follows. 


(1) The caging amplifier for the channel of interest is 
disabled and precession is manually nulled using the 
biasing potentiometer. 

(2) Pickoff and torquer output voltages are recorded. 

(3) A step change is introduced into the torquer through 
the biasing network and simultaneously a stopwatch 
is started. 

(4) The new torquer output voltage is recorded. 

(5) The time required for the pickoff voltage to increase 
to a value much larger than the initial is recorded. 


With the precession rate known, the pickoff scale factor 

is computed as 


K g = AV g /AtAV T K T 

where 


= pickoff scale factor (mV/prad) 

AV = change in pickoff voltage (mV) 

D 

At = elapsed time (sec) 

AV t = change in torquer amplifier output voltage (V) 

K t = torquer output scale factor ((jrad/sec-V) 
o 
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Several measurements of each scale factor were made and 
averaged. The scale factors for the X, Y, Z, and X RD axes are 6.03, 
6.70, 7.51, and 5.74 mV/[irad, respectively. 

(7) Caging Amplifiers . It is necessary to keep the gyro 
floats near null to maintain linear pickoff performance and to prevent 
the floats from drifting into their stops. However, it is also necessary 
to keep caging loop dynamics from interfering with the primary function 
of sensing changes in attitude. To satisfy these needs, rebalance servo 
loops were designed with bandwidths less than the lowest signal frequency 
of interest (0.1 rad/sec). In addition, the original Minuteman caging 
amplifier was kept to control the noncritical redundant axis. 

It was further desired to have the caging amplifiers be 
able to cancel out earth rate without requiring a large position offset. 
For this reason, a proportional plus integral control law was employed 
for the three primary channels. The proportional feedback of the redun- 
dant axis amplifier was not altered. 


(8) Caging Amplifier Gyro Supervision . At low frequen- 
cies, the motion of float relative to the housing can be adequately mod- 
eled as 

• i c d 

(7=0) + to + to 


where 


a = pickoff rate 

go* = inertial angular rate of the housing about the pickoff axis 
to° = commanded rate 
to d = drift rate 


c id 

In steady state, to = -(to + to ). With the pickoff providing a voltage 
V = K 0 , and letting V be the integral of V , the feedback control 

g g I g 

law becomes = -C.V -C V T , with arbitrary gains C and C . This 
F 1 g Z 1 -L £ 

voltage is then put into the torquer amplifier which precesses the float 
at a rate 
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where 


c 


to 


= Vf 


K = overall torquer scale factor 


(The scale factor K is related to the previously defined torquer scale 
factors K t and K>r by K = 92.1 Kt Kt , where 92.1 is the resis- 
tance in ohms of the amplifier current-sensing resistor. ) Combining all 
of the above relations, Laplace transforming, and neglecting initial con- 
ditions results in the equation 

S„ + K T K g C l( r + K^C 2 | = to 1 + io d 


Multiplying through by s or differentiating the original expressions 
yields 

( 2 2\ i rj 

s +2 £( 0 ^ + w)a=s(to + to ) 

where 

w n = i \/ K T K ^2 “ contro1 loop natural frequency 
£ = K^C^/4C 2 = control loop damping ratio 


The gains and C ^ were selected to produce a natural frequency of 

0.1 rad/sec and a damping ratio of 0.7, 

The caging amplifiers were implemented using one opera- 
tional amplifier in each circuit. Figure 6 is a simplified block dia- 
gram of the caging loop, while Fig. 7 shows a basic schematic of the 
caging amplifier. 

An additional small resistor in series with a switch is 
tied between the summing junction and the output of the operational am- 
plifier to allow manual inhibiting of the caging amplifier output. 
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Fig. 6. GST RO CAGING LOOP. When the platform is autocollimator 
stabilized, the gyro pickoff signals can be fed back to the 
torquers to compensate for earth rate and gyro drift. 
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Fig. 7. BASIC CAGING AMPLIFIER MECHANIZATION. A single 
operational amplifier provides both proportional and 
integral control. 


(9) Caging Amplifier Autocollimator Supervision . A second 
caging loop is required when the gyros provide the platform stabilization 
signals. In this case, the gyro pickoff signals are held at null by the 
servo. However, any error torques on the gyro float cause the entire plat 
form to drift away from its original orientation. To maintain constant 
orientation, error signals from the autocollimators are used as inputs to 
the caging amplifiers. The caging amplifiers then cause the torquer am- 
plifier to exactly cancel the drift torques in the steady state. Figure 
8 shows a simplified autocollimator supervision loop. 
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Gimbal Caging 



Fig. 8. AUTOCOLLIMATOR SUPERVISION LOOP. To maintain 
the orientation of the platform while under gyro con- 
trol, autocollimator signals are used to torque the 
gyros. 

(10) Heaters . It is essential to keep the flotation 
fluid at a temperature which makes the float neutrally buoyant in order 
to minimize drift torques. Tight temperature control is maintained by 
a new set of electronics using the temperature sensitive resistor mounted 
on each gyro as the sensor. This resistor is used in a bridge circuit, 
the output of which is compared to a precision reference voltage. An in- 
tegrated circuit thyristor controller is used to switch a thyristor on 
and off, allowing a sufficient number of 60 Hz full and half waves to 
pass through the gyro heater coils to maintain temperature. High gain 
in a small linear region allow the temperature of each gyro to be held 
at 40°C ± 0.1 °c. The heater circuits are completely isolated from the 
rest of the system, and are normally operated continuously to avoid long 
period thermal transients. 
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b. Autocollimators 

Optical reference for the Fixed Base Simulator is provided 
by two different autocollimators: one for maintaining azimuth orientation, 

the other as the primary test sensor. The azimuth unit , which was used in 
a spining vehicle simulator constructed in the Stanford Guidance and Con- 
trol Laboratory [LI], has good wideband performance and moderately low 
noise. Recent advances in solid state electronics, however, made possi- 
ble significant improvements in noise levels and null stability. These 
advances were used advantageously in the design of the primary autocolli- 
mator, and resulted in an improved sensor described below. 

(1) Primary Sensor Configuration. The newly developed 
autocollimator is pictured in Fig. 9. Visible externally are the main 
body, objective lens and mount, light source housing and heat sink, and 
preamplifier box. Also shown is the reference mirror and its adjustable 
mount. 


(2) primary Sensor Optical and Mechanical Arrangement . 

A functional diagram of the autocollimator, Fig. 10, shows the arrange- 
ment of the main mechanical and optical components. 

The infrared emitting diode is switched on and off at 2.2 
KHz. Radiated energy is collected by the field lens and passes through 
a 0.050 in. square aperture located at the focal plane of the objective 
lens. (A square rather than a circular aperture is used to enhance lin- 
earity.) A portion of the emerging energy is diverted 90° by a cubical 
beamsplitter, and passes through the objective lens, emerging as a colli- 
mated beam. The beam is then reflected off the reference mirror, and 
passes back through the objective lens forming a focused image of the 
square aperture on the four quadrant detector as indicated in Fig. 11. 

(3) Primary Sensor Signal Conversion . The position of 

the aperture image on the photodetector is converted into measurements 

of the small angles 6 and 6 between the autocollimator optical axis 

x y 

and the reference mirror normal as discussed below. Each quadrant of the 
diffused barrier photodetector creates a small current proportional to 
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FIG. 9 PRIMARY AUTOCOLLIMATOR AND REFERENCE MIRROR ASSEMBLY. Signals from the 
autocollimator provide two axes of optical reference, A second similar 
autocollimator senses motion about the azimuth (vertical) axis. 



REFERENCE 



Fig. 10. PRIMARY AUTOCOLLIMATOR FUNCTIONAL DIAGRAM. 
Infra-red rays originating at the same point of the 
aperture emerge from the objective lens collimated 
(parallel). After being reflected through an angle 
twice that of the reference mirror, the rays are 
focused on a four-quadrant photodetector. 
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W= 0.0475 in 
9 = 0.005 in 
a = 0.050 In 



ILLUMINATED 

AREA 



PHOTODETECTOR 

SENSITIVE 

AREA 


Fig. 11. GEOMETRY OF THE FOUR-QUADRANT PHOTODETECTOR AND 
APERTURE IMAGE. Displacement of the illuminated area 
from the center of the quadrants results in voltages pro- 
portional to the x and y distances. 

the total energy flux into its surface. In addition, a dark current, 
similar to the saturation current of a transistor, is created indepen- 
dently of the external light. The level of dark current is subject to 
long term drifts. Finally, noise is produced whose power spectral den- 
sity varies inversely with frequency (often called flicker or 1/f noise) 
It is assumed that the energy flux density across the ap- 
erture image is constant, and therefore the total energy flux is propor 
tiorial to the illuminated area. Referring to Fig. 11, the currents pro- 
duced in the four quadrants are 
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‘3 ■ ‘.(t 1 - -X^t 1 - 0 * * "3 

*. ■ ^ * «)(=-h - 0 * \ * ■« 

with |x|, |y | < w - | + | 

where and i^, m = 1,2, 3,4 are the individual dark currents and 

noise, respectively; and is the constant of proportionality, assumed 

to be the same for all quadrants. Each current is converted to a voltage 
in a set of integrated circuit preamplifiers immediately adjacent to the 
photodetector and sent to the main electronics rack. Here, the individ- 
ual voltages are first summed and differenced as follows: 


where 


v x = (V x ♦ v 4 ) - <v a ♦ v 3 > = V ♦ v + n x 

d 


V y = (V x ♦ V 2 ) - (V 3 + V = k y y + V + n y 


V = preamplifier outputs (m = 1,2, 3, 4) 
m 

k = k = proportionality constants 
x y 

V , V = voltages due to dark current 
X d y d 

n^,n^ = noise voltages 


The resulting voltages and are therefore direct 

linear indications of the displacement of the centroid of the aperture 
image. This displacement is further related to angular rotations of the 
mirror relative to the autocollimator by the linearized expressions 


0 = y/2f 

A 

e = -x/2f 
y 
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as can be deduced from the geometry of Figs. 10 and 11. 

The signals V x and V are further amplified and demod- 
ulated. Since the power spectral densities of the dark current and the 
1/f noise are greatest at low frequencies, and the signal occurs at 2.2 
KHz, demodulation effectively eliminates these two error sources. Also 
eliminated are any biases in the amplifiers. The demodulated signal is 
then amplified and actively filtered above 40 Hz to attenuate 2.2 KHz 
ripple and high frequency noise, and to prevent autocollimator mount res- 
onances from destabilizing the platform. The demodulated signal is also 
amplified and filtered above 120 Hz for high bandwidth instrumentation 
purposes. An automatic gain control loop provides additional gain sta- 
bility and linearity by adjusting the output of the IR diode to keep the 
sum of the four quadrant preamplifier outputs at a constant level. 


(4) Primary Sensor Alignment and Coll imat ion . The auto- 
collimator was aligned and collimated using the following procedure. A 
Wild T2 theodolite was first sharply focused at infinity by converting 
it to an autocollimator, pointing it at a mirror, and bringing the re- 
flected image of the projected reticle into focus. The mirror was then 
replaced by the autocollimator and positioned such that the theodolite 
could look into the autocollimator through the objective lens. An ex- 
ternal lamp was placed near the objective lens to illuminate the inter- 
ior of the autocollimator. The detector face was then placed at the fo- 
cal plane of the autocollimator objective lens by adjusting the threaded 
lens mount until a sharply focused image of the detector face appeared 
in the theodolite eyepiece. At this position the image distance of the 
detector face is infinity, enabling the theodolite to focus the detector 
image. 


The IR beam was similarly collimated by bringing the image 
of the square aperture into focus by sliding the aperture plate in and 
out by means of a differential screw assembly. It was necessary to re- 
place the IR emitting diode by a similar diode emitting in the red spec- 
trum in order to make the aperture visible. The degree of collimation 
required made unnecessary a correction for the different wavelength of 
the IR and visible red emitting diodes. 
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Alignment of the detector to the aperture, to center the 
aperture image and make their sides parallel, was accomplished by ob- 
serving the detector and aperture images in the theodolite, and rotat- 
ing and translating the detector in the focal plane with adjusting screws. 
The IR diode was reinstalled and the detector mounting plate, objective 
lens, and aperture plate were locked in place by set screws. 

During the original assembly, the IR diode was defocused 
to ensure uniform illumination of the detector across the aperture image. 
Defocusing was accomplished by adjusting the distance between the diode 
and the field lens so that the image of the diode was focused on the ob- 
jective lens. 

(5) Primary Sensor Performance . Scale factor, linearity, 
noise, and null stability measurements were made using a Davidson D-707 
autocollimator, which has a resolution of 0.1 arcsec. A precision, ad- 
justable, two-sided mirror was placed between the two autocollimators 
such that both were in their linear range. By introducing small changes 
in the orientation of the mirror and recording the outputs of both in- 
struments, scale factor and linearity information was obtained. Noise 
estimates were made using the data recording and processing procedures 
described in Chapter IV. Twenty- four hour null stability was determined 
by continuously displaying the autocollimator outputs on a Sanborn strip 
chart recorder. Results of all these tests are summarized in Table 1. 

Also included are the characteristics of the useful channel of the azi- 
muth autocollimator. 


(6) Azimuth Orientation Sensor . The second autocollima- 
tor, which provides the necessary information to stabilize the platform 
in rotation about the boresight of the precision autocollimator, is also 
a two-axis device; however, one channel is not used. Mechanical arrange- 
ment of the unit is essentially the same as the primary sensor. Basic 
differences of the secondary sensor are: 

(a) incadescent bulb light source 

(b) dc electronics 
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Table 1 


AUTOCOLLIMATOR PERFORMANCE CHARACTERISTICS 



6 

X 

(Control) 

e 

y 

(Control) 

G 

X 

( Instrumentation) 

6 

y 

( Instrumentation) 

G 

z 

(Azimuth) 

Scale Factor 
(mV/ farad ) 

5.70 

5.99 

5.63 

5.92 

7.36 

Linearity Error 

(p 

<0.1 

<0.1 

<0.1 

<0.1 

1 (±2.5 mrad ) 
2.5 ( ±5 mrad ) 

Linear Range 
(mrad ) 

±2.3 

±2.3 

±2.3 

±2.3 

±5 

Standard Deviation 
of Noise Equivalent 
Angle (farad) 

(0.1 to 100 Hz) 

0.451 

0.506 

0.610 

0.550 

2.78 

24 Hour Null 
Stability (farad ) 

<5 

<5 

<5 

<5 

No Data 

Power Required (W) 

4.5 

5.0 














(3) 


four phototransistor detector 


(4) ground glass image focusing surface with four light 
pipes leading to phototransistors 

(5) partially silvered mirror beam splitter 

(6) circular aperture 

Performance of this sensor is shown above in Table 1. 

For more details see [Ll]. 

2. Control Console 

During system layout design, the need for a control center be- 
came evident. As a result, a large mobile 19 in. panel console, pictured 
in Fig. 12 , was employed to provide front panel space and storage. 
Front panels consist of: 

(1) control 

(2) meters — analog and digital 

(3) switchable instrumentation 

(4) hardwired instrumentation 

(5) toroid to platform connectors 

(6) 28 VDC power supply control 

Storage is provided for the following items: 

(1) electronics rack 

(2) 28 VDC power supply 

(3) ±15 VDC power supplies 

(4) gyro spin power isolation transformers 

(5) heater isolation transformers 

(6) relay panel 

(7) wire interconnect panel 
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FIG, 12 CONTROL CONSOLE, All control, switching, and system monitoring 
functions are conveniently handled with this unit. 


a. Electronics Rack 

All of the electronic circuit cards developed for the 
system, with the exceptions of the gyro and autocollimator preamplifi- 
ers, are contained in an electronics rack inside the control console, 
Fig. 13 . Each removable circuit card plugs into a connector. The 
connectors are interwired and also linked to the rest of the system 
through two large wire bundles terminating in 50-pin connectors at the 
interconnect panel. A separate cable connects directly to the external 
autocollimator cable. Table 2 provides a list of the cards and their 
individual functions. 


Table 2 

ELECTRONIC CARD RACK FUNCTIONAL DESCRIPTION 


Card 

Descript ion 

1 

X and Y gyro torquer amplifiers 

2 

Torquer amplifier decoupling networks 

3 

Z and Xpp gyro torquer amplifiers 

4 

Gyro caging (torque rebalance) amplifiers; 
±6.2 VDC precision voltage regulator 

5 . 

Mode switching logic and relay drivers; 
1.32 MHz crystal controlled clock 

6 

Gyro magnet excitation 5 VDC power supply; 
mode switching logic 5 VDC power supply 

7 

5 VDC, 1000 raa autocollimator LED power 
supply; gyro temperature controller 

8 

Primary autocollimator electronics 

9 

Azimuth autocollimator electronics 

10 

Autocollimator instrumentation output 
filters; ±12 VDC power supply 
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FIG. 13 ELECTRONICS RACK. Most of the new electronics required by 
simulator are mounted on cards which plug into this rack. 

The rack is interwired to provide power and signal distribution. 



b. Control Panel 

(1) Arrangement . The simulator is operated from a con- 
trol panel near the top of the control console. As shown in Fig. 14, 
the panel is divided into several areas: 

(1) power switches and display lamps 

(2) mode select switches and display lamps 

(3) caging amplifier inhibit switches 

(4) gyro torquer biasing switches and potentiometers 

(5) autocollimator output biasing switches and potentiometers 

All of these controls except item (2) are self-explanatory or have been 
previously described. 

(2) Mode Switching . The mechanics of switching modes is 
included here; however, discussion of the various modes is deferred until 
the next section. Selection of modes 1 through 6 is accomplished by de- 
pressing the "hold" momentary switch, setting the three mode select 
switches arranged in binary order, and then releasing the "hold" switch. 
Three lamps indicate if the select switches are on. 

Each select switch controls one logic level input to a 
clocked three-bit register. A 1.32 MHz oscillator provides continuous 
clocking unless inhibited either by the selection of disallowed modes 0 
and 7, or by depressing the "hold" switch. A panel alarm light also in- 
dicates if the select switches are in mode 0 or 7. The "hold" switch al- 
lows the operator to set all switches before releasing into the next mode 
and prevents incorrect modes from being entered during switching. The 
oscillator is included so that if a power surge or other anomaly falsely 
switches the register, the register will be reset to the desired mode be- 
fore the mode-switching relays can change state. 

The three bit register output is decoded in a logic cir- 
cuit, which then controls current to the various relay coils through a 
set of transistor relay drivers. The 1.32 MHz oscillator output is also 
counted down to 2.2 KHz in the autocollimator electronics. 
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FIG. 14 


CONTROL PANEL. Except for primary +28Vd-c power switching 
control functions are accomplished on this panel. 


all 



c. System Monitoring 


A desire for convenient and flexible system monitoring 
strongly influenced the configuration of the control console. Figure 
15 displays the arrangement of the system components and cabling and 
emphasizes the advantages of a central monitoring area- 


CONTROL CONSOLE 


r j 



Fig. 15. SYSTEM INTERCONNECT DIAGRAM. in addition to providing 
control and monitoring functions, the control console serves 
as an interface between the Minuteman I electronics and the 
platform. 
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The test point panel consists of 31 connected pairs of 
jacks: one side tied to a signal, the other to ground. The purpose of 

this panel is to provide multiple outputs so that several channels of 
data may be recorded simultaneously. External command signals are also 
introduced to the system through these jacks. 

The monitor switch panel consists of four sets of two ro- 
tary switches with output jacks, enabling any of 61 signals to be moni- 
tored. This panel allows rapid selection of signals to be displayed on 
up to four instruments; e.g., an oscilloscope. All signals appearing at 
the breakout panel are also sent to the monitor switch panel. Channel 4 
of the monitor switch panel is permanently connected to two digital dc 
voltmeters on the meter panel; one with a range of ±2 VDC, the other ±100 
VDC. 

Ten signals were selected for permanent display via analog 
dc voltmeters, also mounted on the meter panel. They are: all four gyro 

pickoff demodulator outputs, all four gyro torque amplifier output volt- 
ages, and two control channel autocollimator outputs, 6 and 6 . These 

x y 

meters are useful indicators of system status during autocollimator acqui- 
sition and gyro spin-up and caging transients. 


3. Minuteman Electronics 

To minimize modification of existing equipment, all Minuteman I 
electronics used were left in the toroid except the pickoff excitation am- 
plifier, which was left in its original location mounted at the base of 
the stable platform. 

Other modules not needed by the simulator were either removed 
or disabled. The computer was completely removed from the system by un- 
plugging several connectors. The computer power supplies were disabled 
by removing the dc to dc converter. Also no longer needed was the power 
discrete switch used by the Minuteman in switching from ground to airborne 
power sources. Level detection, originally used to erect the platform, 
was not needed on the simulator so the level detector converter was re- 
moved from the toroid, and the level detecting bubble levels were not 
mounted on the modified platform. Since the gyro torquing arrangement 
was drastically altered, the gyro tuning network and torque converter 
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were eliminated. The final components removed from the toroid, which 
potentially could be used in future work with the simulator, are the 
digital resolver converters used to indicate the platform gimbal angles. 
To use the resolver converter output, however, would require the develop- 
ment of digital accumulators. Since the resolvers were not immediately 
needed, no effort was made to use them. 

The following toroid modules are used in the simulator. 

(1) Platform Control Amplifiers (3) 

(2) Velocity Meter Control Amplifiers (3) 

(3) Gyro Pickoff Postamplifiers (3) 

(4) Velocity Meter Post amplifiers (3) 

(5) Hybrid Circuit 

(6) Direct Current Filter 

(7) Caging Amplifier 

(8) -28 VDC Power Supply 

(9) 4.8 kHz Power Supply 

(10) 400 Hz Precision Power Supply 

(11) 400 Hz Audio Frequency Amplifier 

(12) Signal Generator 

Power from the 28 VDC supply in the control console is deliv- 
ered to the Minuteman electronics through a special cable which eventu- 
ally terminates at a glass-epoxy board installed at the base of toroid. 
Communication between the toroid and the control console is through four 
cables with connectors at both ends. Three connectors on the toroid, 
originally for the stable platform, already had many of the desired sig- 
nals and were used. A fourth new cable brings out many other signals 
needed by the simulator. Additional coaxial wiring was required inside 
the toroid for this purpose. 
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4. 


Plat form Assembly 


The primary purpose of the platform assembly is to provide 
controlled motion for the sensors. The platform assembly, consisting 
of the stable platform, rotary table, and mirror mount, has already been 
pictured in Fig. 2 at the beginning of this chapter. The rotary table 
provides flexible and convenient orientation of the platform control axes 
relative to the directions of gravity and the earth's spin axis. Optical 
reference is provided by the mirror mount, which is rigidly attached to 
the rotary table head. 

a. Stable Platform 

The "inside-out” gimbal arrangement of the Minuteman I 
stable platform proved to be quite amenable to modifications. The gimbal 
arrangement and axes convention are shown in Fig. 16. All three gimbal 


INTERMEDIATE 

GIMBAL 



Fig. 16. PLATOFRM GIMBAL CONFIGURATION. The 
orientation sensors are clamped rigidly to the 
inner gimbal shaft. 
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shafts are supported by ball bearings. Torques between gimbals are ap- 
plied by permanent magnet dc motors whose static gains are 0.027 N-m/V. 
The motor input voltages are provided by the power amplifier stages of 
the platform control amplifiers. The gains of these power amplifiers 
are 302, 164, and 333 for the X, Y, and Z axes, respectively. Command 
voltage inputs to the power amplifier may be applied through dropping 
resistors, giving gains of 2.26, 1.23, and 2.49. 

The original gyro and velocity meter mounting brackets 
were replaced by new mounts for the Atlas gyros and autocollimator. The 
alignment unit, consisting of bubble levels and mirror mounted on a quartz 
block, and supporting electronic components were removed with the original 
brackets. The two gyro brackets are rigidly clamped to the inner gimbal 
shaft and are tied together by the autocollimator and connector mounting 
bracket as pictured in Fig. 17. Several additional brackets were fab- 
ricated to mount the azimuth reference mirror, primary autocollimator, 
balancing weights, and isolation transformers. Three digital gimbal 
angle resolvers were left in place on the platform, although not used 
because of their inferior small angle accuracy. 

Mechanical stops prevent the outer and intermediate gim- 
bals (Z and X) from rotating more than ±110° and ±25°, respectively. 
Inner gimbal travel must be held to within ±90° to avoid pinching wire 
cables. 


b. Rotary Table 

The platform base is rigidly attached to a 24 in. tilting 
rotary table by means of a special adaptor. This hollow, cylindrical fix- 
ture allows passage of external cabling into the platforms and also pro- 
vides room for the pickoff excitation amplifier suspended from the plat- 
form base. 

Elevation angles about an East-West axis between 0 and 90° 
can be set by cranking the tilting handwheel. A vernier scale allows 
reading of the elevation angle to 1 arcmin. The platform may also be 
rotated by crank through any angle about the spindle axis perpendicular 
to the rotary table face. This angle can be repeatedly read out to an 
accuracy of 2 arcsec. With these two-degrees-of-freedom, the outer (Z) 
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AUTOCOLLIMATOR AND 
CONNECTOR BRACKET 
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CO 



GYRO 

PREAMPLIFIER 


X-Y GYRO 
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Y GIMBAL ANGLE 
RESOLVER 


Y GIMBAL 
MOTOR 


X-Z GYRO 
BRACKET 


X GIMBAL 
ANGLE RESOLVER 


FIG. 17 SENSOR MOUNTING CONFIGURATION. New gyro brackets are mated to the 
original cross (autocollimator) bracket and are clamped to the 
inner (Y) gimbal shaft. This arrangement provides rigidity and 
reduces deflections due to the cantilevered gyros. 


gimbal axis can be oriented quickly to any position between vertical and 
horizontal, including earth polar; and the sensors then rotated about 
this axis. 

The table base rests on three adjustable bolts so that 
the table can be leveled. There are also four brackets holding bolts 
located at the bottom corners of the base. These bolts push against the 
base and may be adjusted in various combinations to rotate the table in 
azimuth. Tightening all bolts locks the table base to the floor plate. 


c. Mirror Mount 

A special mount was built to hold the primary autocolli- 
mator reference mirror and the azimuth autocollimator. The crescent 
shaped arm, which is pivoted at the midpoint, curves up and around the 
stable platform to hold the reference mirror directly above the primary 
autocollimator. Although normally locked in the unpivoted position, the 
arm can also be rotated and locked at angles of 10° and 57°. These posi- 
tions allow the platform to be operated with the gimbals nonorthogonal , 
while still maintaining an optical reference. 

\ 

The azimuth autocollimator is installed in a hole in the 
arm below the pivot point thereby remaining permanently fixed relative 
to the rotary table head. This autocollimator reflects off a mirror 
mounted on a compensating weight bracket attached to one end of the in- 
ner (Y) gimbal structure. It therefore senses platform motion about 
the intermediate and outer (X and Z) gimbal axes, although the X 
axis sensing is redundant and is not used. 

The primary mirror mount is cantilevered from the upper 
end of the arm. The mirror is suspended in a two-axis gimbal and is 
adjusted by two micrometers. 


d . Alignment 

Alignment of the gyros and autocollimators involved sev- 
eral stages of adjustments. The sequence of relative alignments used is 
as follows. 
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The rotary table was leveled. 


The two sided primary mirror was made parallel to the 
rotary table head. A Wild T2 theodolite, used as an 
autocollimator, was pointed at the back side of the 
primary mirror. The table head was elevated about 60° 
for this purpose. The table head (and thereby the mirror) 
was then rotated and the motion of the reflected reticle 
image observed. (The origin of the reticle image de- 
scribes a circle if the mirror is not perpendicular to 
the rotation axis.) The micrometers were then adjusted 
to minimize the diameter of the described circle. The 
procedure was iterated until no apparent motion of the 
reticle image was observed. 

The rotary table tilt shaft was oriented East-West by 
referring to the star Polaris, again using the theodo- 
lite. The rotary table, which had already been oriented 
approximately East-West under a door through which Pol- 
aris could be observed, was elevated to the local lati- 
tude. The theodolite was then positioned along the line 
from the primary reference mirror to Polaris which was 
precisely centered in the theodolite. Corrections to 
true North from Polaris were computed from American Eph - 
emeris and Nautical Almanac Pole Star tables, and the 
theodolite moved to the true Polar position. The theo- 
dolite was then plunged 180° about its elevation axis, 
which brought it to an orientation looking at the approx- 
imately polar oriented mirror. The theodolite was con- 
verted to an autocollimator, and the vertical lines of 
the actual and the image reticle were brought into coin- 
cidence by adjusting the table base plate bolts. The 
horizontal lines were brought together by cranking the 
table elevation control wheel. At this point the refer- 
ence mirror was very precisely (error <50 |jrad) perpen- 
dicular to the earth's spin axis and the elevation shaft 
horizontal East -West. 

The primary autocollimator optical axis was made parallel 
to the X-Y gyro spin reference axis. The platform was 
oscillated about the Z axis while under autocollimator 
stabilization control by driving the Z gimbal motor 
power amplifier with a 1 Hz sine wave from a low frequency 
function generator. Since the X and Y autocollimator 
outputs were held at null by the servo loop, sinusoidal 
outputs from the X-Y gyro indicated that the spin axis 
was not parallel to the autocollimator optical axis. 

Screws on the autocollimator mount were then adjusted to 
minimize the gyro outputs. 
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Attempts to align precisely the X-Y gyro and autocolli- 
mator about the Z axis were thwarted by the viscosity of the gyro flo- 
tation fluid. A rotation about the X axis, for example, is accompanied 
by a rotation of the float about the Y axis. This is easily shown by 
the following simplified analysis: 



where 

c = viscous friction coefficient 
x 

T = viscous friction torque 
x 

h = gyro angular momentum 

a , a = float angular rates 
x y 


Values for h/c are presented in Section C. 1. a. 

The gyro outputs are clearly coupled. Since careful at- 
tention had been paid to the original mechanical alignment, it was felt 
that the existing accuracy was adequate. By the same reasoning, no at- 
tempt was made to improve the alignment of the noncritical X-Z gyro 
and the azimuth autocollimator. 

In summary, orienting the reference mirror and azimuth 
autocollimator using the rotary table controls provides the basic labo- 
ratory and inertial reference. The servo loops then hold the sensor 
axes to this reference. 


5 . Power 

Eleven power supplies, four regulators, two crystal oscilla- 
tors, and a motor generator supply power to the system. These in turn 
receive their power from both single-phase and three-phase laboratory 
60 Hz lines. Power distribution is indicated in Figs. 18a and 18b, 
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115 V 60 Hz Line 



Fig. 18b. CONTROL CONSOLE AND PLATFORM POWER DISTRIBUTION. 






























As can be seen from these two figures, the original Minuteman 
toroid and the new electronics have essentially separate power distribu- 
tion systems. They are tied to each other only through the ground buss 
in the control console and through the relay drivers which control the 
flow of 28 VDC power through the relay coils. 


6. Data Acquisition 

a. Strip Chart Recording 

Permanent records of both analog and digital data are 
made using a variety of equipment. An eight channel Sanborn strip chart 
recorder is used to display analog signals directly. This instrument is 
primarily used in obtaining qualitative information, since it is inher- 
ently limited to about 1 percent accuracy by the resolution of grid paper. 
More precise data in a form suitable for digital processing is obtained 
by a procedure illustrated in Fig. 19. 



Fig. 19. DATA ACQUISITION PROCEDURE. Three stages of data 
acquisition and editing are required to put data in a form 
suitable for digital processing. 
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b. FM Analog Recording 


Up to eight signals of interest are tapped at the control 
console test point panel and sent to the strip chart recorder for ampli- 
fication and display. The amplified signals are then recorded on magnetic 
tape by a high fidelity Ampex FR-1300 FM recorder/ reproducer . The amount 
of amplification is selectable and is usually chosen so that the maximum 
output voltages are approximately ±1 V. Amplification is necessary to 
maximize the signal-to-noise ratio in the recorder, which has a maximum 
input level of ±2 V. Since the recorded signals are to be digitized, it 
is necessary to simultaneously record a sine wave from an oscillator set 
at the desired sampling frequency. Recording length is determined by the 
desired number of sampled data points and the sampling frequency. 

c. Digitization and Storage 

After the test run is completed, the recorder and tape are 
moved from the laboratory to an XDS Sigma 5 computer facility located in 
the same building. The appropriate signal playback channels are connected 
to three-pole active low-pass filters to attenuate signal components at 
frequencies greater than approximately 30 percent of the sampling fre- 
quency. This is necessary to prevent high frequencies being aliased into 
the low frequency spectrum. The analog data is normally recorded at a 
tape speed of 3 3/4 ips, but is played back at either 15 or 30 ips, re- 
sulting in digitization time savings of factors of 4 or 8. The increased 
speed requires, however, that the anti-aliasing filter cutoff frequency 
be increased by the same factor. 

The oscillator channel is connected directly to a pulse 
shaping circuit that generates a 1 u sec square pulse at every positive 
zero crossing of the input sine wave. This pulse is then used to trigger 
a multiplexed input analog-to-digital converter (ADC). The advantage of 
triggering from a frequency signal recorded simultaneously with the data 
signals is the elimination of frequency distortions caused by variations 
in tape speed during playback. 

Upon receipt of the triggering pulse, the input channels 
are sampled and digitized at a 140 kHz repetition rate, resulting in a 
14 bit plus sign word for each sampled channel to be stored sequentially 
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in a double buffer. Buffer length is set such that n channels of 256 
words per channel fill each buffer. The number of input channels n is 
variable. Once one buffer is filled, the other buffer is used. While 
the second buffer is being filled, the first buffer is written as a rec- 
ord onto magnetic tape. The buffer roles are then alternately swapped 
until a preset number of digital tape records have been written. 

The format of the digital tape is not convenient for pro- 
cessing since channel data are interlaced. To make the tape format more 
compatible, an editing program is executed which passes through the data 
record by record, selecting points from one channel only, and writing 
out 256 word records until the desired number of records have been writ- 
ten. A file mark is then placed on the tape. This is repeated for a 
preselected number of original channels. 

D. Operation 

The simulator can be operated in any of six modes. As discussed 
previously, the modes are selected on the console control panel through 
twelve triple-pole, double-throw relays. A functional description of 
each of the modes is given below. Diagrams of the relay connections are 
contained in Appendix B. 

Mode 0: Not allowed . 

Mode 1 : Gyro Off 

a. Autocollimator stabilization 

b. No gyro spin power 

c. Caging amplifiers inhibited 

d. Gyro torquer biasing 

Mode 2 : Gyro Spin-up 

a. Autocolliraator stabilization 

b. Gyro spin power 

c. Caging amplifiers inhibited 

d. Gyro torquer biasing 

Mode 3: Autocollimator Control 

a. Autocollimator Stabilization 

b. Gyro spin power 
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c. Gyro caging via gyro pickoff 

d. Gyro torquer biasing 

Mode 4: Gyro Control with 4-Channel Caging 


a. 

Gyro 

stabilization 

b. 

Gyro 

spin power 

c. 

Gyro 

caging via autocollimator (4 axes) 

d. 

Gyro 

torquer biasing 


Mode 5: Gyro Control with 3-Channel Caging 

a. Gyro stabilization 

b. Gyro spin power 

c. Gyro caging via autocollimator (no z channel) 

d. Gyro torquer biasing 

Mode 6: External Stabilization 

a. Stabilization by external compensation 

b. Gyro spin power 

c. Gyro caging via external signals 

d. No gyro torquer biasing 


The simulator can be started in modes 1, 2, 3, or 6. The gyro con- 
trolled modes 4 and 5 require that the gyros be at operating speed and 
caged before entering. Mode 5 is primarily used to keep the platform 
inertially fixed in the polar orientation while maintaining optical refer- 
ence in the polar direction. Mode 6 is to be used when external compen- 
sation is available from a computer or other network which combines the 
gyro and autocollimator error signals. 
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Chapter III 


SIMULATOR MODEL 


A mathematical model of the Fixed Base Simulator is now developed. 
Major system components are first analysed and then combined to evaluate 
the responses of the autocollimator and gyro stabilization loops. This 
model is used in Chapter V to help interpret experimental data obtained 
from the orientation sensors and to develop refined models of the system. 


A. Gyro Model 


The two-degree-of-f reedom Atlas missile gyros can be represented 
as in Fig. 20. The two flywheels of the gyro are treated as one combined 
wheel. 



Fig. 20. FUNCTIONAL ARRANGEMENT OF ATLAS MISSILE TWO- 
DEGREE- OF- FREEDOM GSTRO. The spherical float and the 
gimbal are floated inside the case. Support wires 
center the float. 
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The float and gimbal are restricted to small angles about their null 


\ 


position and are oriented as in Fig. 21a and b. 



<3o 


Fig. 21a. GYRO FLOAT SMALL 
ANGLE ORIENTATION. 


Fig. 21b. GYRO GIMBAL SMALL 
ANGLE ORIENTATION. 


The subscripts i, o, s refer to the inner, outer and spin axes respec- 
tively, and the subscripts c, f, g denote case, float, and gimbal. 

The total angular momentum of the float is 


H. = I 1 + h (3. 1) 

f f w 

where 1^ is the combined moment of inertia of the float and wheel, and 
h is the angular momentum of the wheel with respect to the float. (Nota- 
tion is defined in the List of Symbols, p. xiv. ) Taking the time derivative 
of (3. 1) in the inertial reference frame yields 

i - i • 

H f = IfW " + w _1 X H f , (3.2) 


where it is assumed that the wheel speed is constant; i.e. 

f_ 

h =0 
w 


(3.3) 
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h»mIh- 


Since the float angles are small, 


— f-i — — c- i 

OJ = 0 + W 


(3.4) 


where 


a = 


o 


a. 

1 


Expanding (3.2) using (3.4) 


- fi f Vi 


— c-f — c- 


H f = I fL° + W + ^ x 


. r i. • _ - 

+ (o + a) C 1 ) x I- (a + w C X ) + h 
„ X w 


(3.5) 

Equation (3.5) may now be linearized by assuming that products of angular 
velocities are negligibly small. 


i = ff £ _ = c 

H = 1-0 + 0 X h + F /- 1 + C 0 C x h 
f f w f w 


(3.6) 


Coordinatizing (3.6) in the float reference frame yields 


1-8 +c. h 
f t o i 


I. 5 - o h 
f t 1 ° 


*f a s 
s 


- , • c-i .c-i . c-i ~ 

\ ( “o +<J s W i 'Vs > 

V 


1 (^' i -o s ^“ i + aA C " i > 

1 So o b 


•c-i .c-i >c-i. 
I- ( w «5 + c. w - CJ W . ) 

f s 1 o o 1 

l_ s 


<w®“ x - - (s JO ®' 1 + a w s’ 1)h 
1 so o ® 


, c-i c-i c-i , 
“ (W o +0 s W i " a i% )h 


(3.7) 


where 


^ = combined wheel and float transverse moment of inertia 

t 

p = combined wheel and float spin axis moment of inertia 
s 

h' = magnitude of wheel angular momentum 
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and 


U) 


c -1 


C— 1 


w 


C-1 


GJ 


C-1 


Considering all products of float angles and case angular velocities 
and accelerations as negligibly small, (3.7) may be simplified as 


i 



... • c-i. , c-jL . 

L (a + w ) + (a. + w. ) h 

f O O IX 

t 

(&'i + O ' (i o + <C> h 

t 

... . c-i x 

I (ct + w ) 
vu s s ' 
s 


(3.8) 


With the small angle assumption, the spin axis equation is uncoupled 
from the two transverse axes equations. 

In a similar manner, the gimbal angular momentum is 


H = I w g_1 (3.9) 

g g 


Again, taking the time derivative in the inertial reference frame, 


H = I w g 1 + G) g ” 1 x I w g-1 (3. 10) 

g g g 

Assuming the gimbal angles to be small allows the gimbal angular velocity 
to be written as 


= e 




g 


(3.11) 


where 
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(3. 16b) 


M f = - VV^ > " k i ( V 6 g } + M f 

i i i fe i i 


M = -D 0 +D.ct - T, (0 -a ) - k 0 + M (3.17a) 

g o i o i o o g 


M = -D 0 + D.a. - T 0 - k.(0 -a.) + M t 

g. o g. li o g. l g l g (3. 17b) 


c c 

Mf Q , = torque generator command moments 

D , D = flotation fluid viscous friction coefficients 
o 1 

k QJ k ^ ~ support wire rotational elastic restraint coefficients 

T 0 , T ± = support wire tension coefficients 

m 3 , = float disturbance moments 

1 o x i 

M d , M d = gimbal disturbance moments 
g o g i 

Expanding (3.14) and (3.15), and Laplace transforming results in 


l f s +D 1 s+T i 

t 


-(D i S+T i ) 


I f S +D i S+k 1 


-(D i S+k i > a i <S) 


-(D i S+T i ) 


I S +D S+ (T.+k ) 
g t ° 1 ° 


e_ cs) 


-<D i s+k i ) 


I S +D S+(T +k ) 0 (S)j 

g t O o i g t 


I f S h 
t 

~h If t S 

IS 0 
s t 


C 1(s > 




(S) 

o 

* Y (s) 

“k (S) 
8 o 

M? <S) 


(3.18) 
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Initial conditions are neglected since only steady state response is 
desired. 

Parameter values for the gyros are listed in Table 3. 


Table 3 

EXPERIMENTALLY EVALUATED GYRO PARAMETERS 


Parameter 

X-Y Gyro 

X-Z Gyro 

Units 

h 

1.0 

1.0 

2 

Kg-m /sec 

D i 

0.048 

0.028 

N-m-sec 

D 

0. 51 

0.41 

N-m-sec 

o 



k. 

l 

0.0014 

0.00093 

N-m/rad 

k 

o 

0.0024 

0.0019 

N-m/rad 


0.0016 

0.0016 

„ 2 
Kg-m 

I 

0.022 

0.022 

2 

Kg-m 

^t 



T i 

1.35 

1.35 

N-m/rad 

T 

1.35 

1.35 

N-m/rad 

o 




The design value for the angular momentum h is used since the wheel 

moment of inertia is unknown due to lack of documentation. Values for 

the viscous dampling coefficients , D q and support wire rotational 

elastic restraint coefficients k., k and the tests used to obtain them 

1 o 

are discussed in Chapter II. These values are accurate to within 10 $. 

The float and gimbal inertias and the support wiz'e tension coefficients 
were derived from tests performed on the X-Y gyro with the spin power 
off. Sinusoidal torques at numerous frequencies were applied to each 
axis of the float and the pickoff voltages recorded. The transfer func- 
tions between the input torques and pickoff signals with zero angular 
momentum were derived and fit to the experimental frequency response. 

Large uncertainties, perhaps as high as 100$, exist in the inertia and 
tension estimates due to pickoff noise and variability of the experimental 
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data. However, it will be shown later that these parameters have an 
insignificant effect on overall system response, and as a result no 
effort was made to improve the estimates. Values of the inertias and 
support wire tension coefficients for the X-Z gyro are assumed to be 
equal to those of the X-Y gyro, since the gyros are of identical construc- 
tion. These values also are of little consequence in the ensuing dis- 
cussion. 

Several additional approximations and assumptions will now be made 
in order to simplify (3.18). 

1. k Q , are negligibly small. 

2. D » D. 

o 1 

3. The inner gimbal equation is decoupled from the rest of 
the system due to smallness of the term (D^s + k.) 

It will no longer be considered. 

4. The case perturbations are limited to small angles 6 C 

and 0 C± . ° 

5. The command torques, derived from the caging ampli- 

fiers exactly cancel the constant earth rate term and any 
constant terms in the disturbance torques. 

6. Caging amplifier dynamics are sufficiently low in fre- 
quency to be neglected. 

7. The angular velocity of the case u) C-i consists of a 
constant term due to rotation of the earth and a zero 
mean perturbation due to motion of the gyro mount rela- 
tive to the earth. 


With these assumptions (3.18) becomes 



(3.19) 


where and Mf are the residual float disturbance torques after 

the mean is removed. 
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The roots of the characteristic equation associated with (3. 19) are 
well separated, and, as a result, the characteristic equation can be factored 
approximately into 

°' !0> 



where 



(625 


rad/sec) 


2 



0.07 


T. 

u) = — - =2.6 rad/sec 
o 

D 

co = — — = 23 rad/sec 
2 I 

g t 

The approximate transfer functions between the case angles 0 and the 

pickoff angles o are presented in Table 4. 

From these transfer functions it can be seen that the only major 

departure of the Atlas gyro from direct angle indication is the large 

low frequency coupling of outer axis case motion into inner axis pickoff 

angle (o./0 o = 0.51, w < 1) . This coupling was acceptable in the orig- 

c 

inal Atlas platform since the pickoff angles were held at null by the 
platform servo loops. 


B. Autocollimator 

As discussed in Chapter II, the primary autocollimator has two out- 
puts, one for control and the other for instrumentation. The transfer 
function between base angle and control channel signal is 



K. 


AC 


V + JT- 

oj . _ AC 

AC 


+ 1 


(3.21) 
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Table 4 


RESPONSE OF GYRO FLOAT AND GIMBAL TO CASE MOTION 


Translcr 

Function 

Paramol rlc 

Numerical (X-Y pvro) 

Remarks 

C o ' 
O 

C 

. w? s [’'/ * (^r * § ‘ 1 


S 1 



k 

\ w n n / 

0.040 

S 2 2 (0.07>S . , 

5 + ~bIs — ♦ 1 

(625) 

snail 

a t 

5“ 

°c 

n (m; s * ■) 


large cross-coupling 
1 S| « 3 

IB 

, V* ft s * l ) ! 

0.001®(^ * i) 2 

S 1 

■'(* ■?•■)£• $«•■) 


e 

*o 

o 

c 



&»* 



e 

g o 

V 

C 


°' o48 (A * i > 

snail 



m 

0 

0 

negligible 

c 

HbSHH 

s a 

-0.016 — - =— 

(o *•)(£♦■> 

»0 | s| « 2.0 
« 1 |0| > 23 

































where 


7 AC’ 7 AC 

x y 
e , e 

x’ y 


K AC x ' K ACy 


co 


AC 


= autocollimator control output signals (V) 
= autocollimator base angles (rad) 

= 5700, 5990 V/rad 
= 264 rad/sec 


f; = 0.71 


The transfer function for the instrumentation channel is 


where 


AI 

e 



(3.22) 


AI x' 


l AI 1 


X 


7 = autocollimator instrumentation output signals (V) 

AJ. 

y 

K = 5630, 5920 V/rad 

AJ. 

y 

cu =755 rad/sec 

All 

u) =1250 rad/sec 

ai 2 

£ =0.71 


Due to flexibility in the autocollimator mount, the base angle is 
composed of both the " rigid” displacement of the platform gimbals and 
the bending angle of the mount. The latter is considered a disturbance 
input to the autocollimator. 


C. Platform 

The Fixed Base Simulator is normally operated with the gimbals 
orthogonal, although it can be stabilized under gyro control with the 
gimbals non-orthogonal. Since all data were taken with the gimbals 
orthogonal, the platform model will be limited to this case. Deriva- 
tions of the general non-orthogonal three-axis platform equations can 
be found in [W7] and [B9]. 
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The linearized rotational equations of motion of the platform with 


orthogonal gimbals are 


: e = m + 

P P rn 


M. + 

D 


M + 
gr 


M 


(3.23) 


whe re 


M 


m 


M 


b 


M 


gr 

M , 


= angular acceleration of gimbals relative to base 
= effective moment of inertia matrix 
= gimbal motor torque 
= gimbal bearing torque 
= gyro reaction torque 
= disturbance torque 


Values of the on-diagonal moments of inertia obtained from a closed 

loop frequency response of autocollimator output to command torque input 

are: I =0.70 Kg-m 2 , I =0.25 Kg-m 2 , I =0.52 Kg-m 2 . The effects 
xx yy zz 

of off-diagonal terms are small and are treated as disturbances. 


D. Motor Torque 

The torque from each of the three permanent magnet gimbal motors is 
proportional to the armature current, 

M = K i (3.24) 

mm v 7 


where 

K m = motor torque gain (0.53 N-ra/A) 
i = armature current 

The current, in turn, is related to the platform control amplifier out- 
put voltage and the back EMF of the motor by 
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i 


(3.25) 



whe re 

R = armature resistance (23 fi) 
a 

L = armature inductance (23 mH) 
a 

V = control amplifier output voltage 
P 

E = back EMF 

u> = R /L = 1000 rad/sec 
m a a 


The motor torque gain was deduced from measurements of the torque re- 
quired to balance motor torque at several values of input current. The 
torque was supplied by a calibrated spring balance acting through a known 
moment arm. Errors introduced from uncertainty in the moments arm and 
calibration of the spring balance were limited to less than 5 e f. 

The back EMF of the motor is 


E = K b 0 p (3.26) 

where 

K = back EMF proportionality factor (0.53 V-sec) 
b 

In order to satisfy the condition that input power must equal output 

power plus power lost, the back EMF coefficient must equal the motor 

torque gain. An independent check of the back EMF coefficient was made 

by oscillating the platform by hand and recording the voltage induced 

across the motor armature. Inis low accuracy experiment produced a 

value of K. = 0.6, well within the limits of experimental error, 
b 

The output of each platform control amplifier is 


V = K (V 
P pc 


K d V T 


(3.27) 
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where 


V 

c 


= output of compensation network 



= command torque voltage 



= platform control power amplifier gain 
= command torque gain 


The values of the X, Y, and Z axes amplifier gain 
V Dc/ V DC' three comjnand torque gains are 0.0075 V DC /V DC 

these electronic measurements is approximately Vh. 


K p are 302,164,333 


Accuracy of 


E. Compensation Network 

Each Minuteman I platform control amplifier has a passive compensa- 
tion network. Its input is either an autocollimator output or a demodu- 
lated gyro pickoff signal. The output is amplified by a power stage to 
drive the gimbal motors. 

The compensation can be separated into three main functions: lag- 

lead to increase low frequency gain, lead-lag to provide damping and 
stability margins, and high frequency noise attenuation. The transfer 
function, computed from a network schematic, is 



where the break frequencies in rad/sec are 


w - 0. 15 w =30 w = 1200 
C 1 C 3 C 5 


co =10 
C 2 


w = 200 w = 5000 
C 4 C 6 


Fig. 22 is a Bode plot of the compensation network. 
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F. Gyro Reaction Torque 

Motion of the platform gimbals is resisted by the two gyros through 
tension in the gimbal support wires and viscous coupling between the gyro 
gimbals and case. From (3. 17), throwing out unnecessary or small terms, 
the torque applied by the gimbal to the case is 


M = D e (3.29) 

Sr c ° g G 


M = D 6 +T0 

g o g. o g. 

fa r i i & i 


(3.30) 


The two gyros are oriented such that both the outer gimbal axes are 
directed along the platform X axis, Fig. 23. The inner gimbal axes 
of the x-Y and X-Z gyros are in the direction of the platform Y and 
Z axes, respectively. 
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(3. 33) 



D 0 + T 9 

o g. o g 

2 2 


where the subscripts 1 and 2 refer to the X-Y and X-Z gyros, respec- 
tively. 


G. Gimbal Bearing Torque 

The gimbal shafts are supported by ball bearings which provide close 
alignment and low friction between the shafts and their supports. The 
friction normally is considered to be caused by a combination of rolling 
and sliding resistance, and the level of friction to be independent of 
the relative velocity and displacement. However, for the small motion 
exhibited by the simulator it was found that a more descriptive model 
is a lightly damped rotational spring. The background and supporting 
evidence for this result is presented in Chapter V. It is an important 
example of how power spectral density results can be used in developing 
system models. 

Nonlinearity of the friction in the restraint torque makes analysis 
of the system difficult. However, for small motion of the platform gim- 
bals, the bearing coulomb friction is small compared to the viscous 
damping provided by the gyros and can be neglected. The gimbal bearing 
torque is therefore modeled as 

M, = -K 9 (3.34) 

b s p 

Experimentally determined values of the X and Y axes elastic 

restraint coefficients are K_ =410 N-m/rad and K_ = 380 N-m/rad. 

x s y 

The Z axis coefficient was not measured since only the primary X and 
Y axes are considered in the experiments of Chapter V. 


H. Closed Loop Dynamics 

The equations describing the major system components are now com- 
bined so that the closed loop behavior of the system can be described. 
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The simulator can be stabilized using either the gyro or autocollimator 
error signals, and these two cases must be treated separately. 


1, Autocollimator Stabilization 

Using the autocollimators as the reference, the simulator stabi- 
lization loops separate into three independent channels. Products of 
inertia and the elastic restraint terms of the gyro reaction torques are 
small and are considered disturbances. Figure 24 shows the functional 
relationships of the major components. 

The figure as shown represents only the Y platform axis. In 
the X axis, both gyro gimbals provide damping to the platform, and there- 
fore a second parallel gyro gimbal feedback path around the platform block 
should be included. The transfer function of these blocks must also be 

changed to 6 /Q 0 (Table 4). Gyro pickoff cross-coupling indicated 

b o c 

at the bottom of the figure can be neglected for the X axis. 

The Z axis block diagram is similar to the Y axis shown, 
except that the Z autocollimator does not have an instrumentation channel 
or noise filtering. The autocolliraator block here is more appropriately 
modeled as a constant gain K . 

In all three channels, motor back EMF can be neglected. The motor 
lag and the two highest frequency lags in the compensation network can also 
be neglected when considering frequencies less than 15 Hz. Above this fre- 
quency, phase shifts from these elements become important. 

Further noted in Fig. 24 are the sources of disturbance and com- 
mand to the system. Starting at the upper left of the diagram, the auto- 
collimator output can be biased by introducing a signal 7 . Biased and 

c 

unbiased autocollimator outputs are then physically available at the points 


denoted 7 ^ and 7 ^, respectively. Platform torque commands ai'e intro- 
duced by applying a voltage Vm to a divider circuit with gain K. at 

A d d 

the input to the platform motor power amplifier. Environmental and other 

unspecified torques directly accelerate the platform gimbals. All 

gyro gimbal disturbances are included in the small additive term n . 

G 

Autocollimator noise n^ is created primarily in the photodetector and 
in the flexible autocollimator mount. This noise is filtered in both the 
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Gyro 

Gimbal 



Fig. 24. SINGLE-AXIS AUTOCOLLIMATOR STABILIZATION LOOP. At mid-frequenices (0.1 to 100 Hz), 
back EMF, motor lag, and gyro gimbal dynamics have a negligible effect on the response 
of the loop. 












control and instrumentation output channels of the primary X-Y autocol- 
limator. Finally, gyro electronic noise and the effects of disturbance 


torques are combined in the additive input n^. 

It will be of later interest in evaluating experimental results 
to have expressions for the closed loop transfer functions between the 
sensor outputs and the various command and disturbance inputs. As a first 
step, a transfer function for the platform section is derived. 

The transfer function between the motor torque M and the 

m 

platform angle 0^ including the gyro gimbal dynamics is almost identi- 
cal to the transfer function obtained by replacing the gyro dynamics with 
a gain of unity. In other words, the gyro gimbal dynamics are not im- 
portant. The reduced transfer function for the Y axis is then simply 



I S + D S + K 

yy o s 


(3.35) 


The X axis damping coefficient is twice as large since both gyros are 
viscously coupled. 


x ~ 


M 


m 


I S + 2D S + K 
xx os 


x 


(3.36) 


Figure 25a is a root locus of the autocollimator stabilization 
loops as a function of autocollimator gain. The actual operating gain 
is sufficiently low that the poles of the closed loop are quite close to 
those of the open loop. 

An enlarged view of the locus near the origin is shown for the 
X and Y axes in Fig. 25b and 25c. Note that the presence of the auto- 
collimator filter poles at (-180, ± jl80) sec * prevents the roots orig- 
inating at the lower frequency platform gimbal/bearing poles [X = (-0.8, 
- j24.2) sec and Y = (-1.1, ± j39.0) sec *] from penetrating very 
far into the left half plane. As a consequence, the closed loop poles 
lying on this branch of the locus are lightly damped. 
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FIG. 25a AUTOCOLLIMATOR STABILIZATION LOOP OVERALL ROOT LOCUS 



FIB. 25b X-AXIS AUTOCOLLIMATOR STABILIZATION LOOP LOW FREQUENCY ROOT LOCUS 



I m( S) 



rad /sec 


FIG. 25c Y-AXIS AUTOCOLLIMATOR STABILIZATION LOOP LOW FREQUENCY ROOT 
LOCUS 


y AXIS LOCUS 





Closed loop transfer functions between most of the system 
inputs and outputs indicated in Fig. 24 are expressed in Table 5 . Roots 
above 100 Hz are of little interest and are neglected for simplicity. The 
positions of the closed loop poles are indicated in the root loci of Fig. 
25b and c. 

The transfer functions between autocollimator noise n and 
the sensor outputs will be of particular interest in the discussion of 
experimental results in Chapter V. Magnitude plots of these transfer 
functions, based on the parameter estimates indicated in this chapter, 
are shown in Fig. 26. 

2. Gyro Stabilization 

Using the gyros as the reference, the simulator stabilization 

loops are coupled due to the gyro viscosity and cannot be treated as 

three independent channels except at higher frequencies. The coupling, 

however, is almost entirely from the X axis into the Y and Z axis 

and not in the reverse direction. Figure 27 is a block diagram of the 

coupled X and Y stabilization loops. The unidirectional coupling 

is indicated by the block a. (S)/0„ (S) evaluated in Table 4 . The mag- 

i °c 

nitude of this transfer function is less than 0.1 above 20 r/sec falling 
off inversely with increasing frequency. 

At higher frequencies, where coupling can be disregarded, the 
individual stabilization loops are nearly identical to the autocollimator 
stabilization loops except that the gyro does not have filtering. Fil- 
tering is not required in the gyros since their mounting brackets are 
sufficiently rigid to prevent instabilities. 

The closed loop poles of the X-Y system in Fig. 27 factor, by 
inspection, into the poles of each individual axis and the poles of the 
coupling transfer function. Root loci can therefore be drawn for each 
axis; Fig. 28 shows the loci for the X and Y axes. The gyro loops 
are characterized by improved damping of the dominant oscillatory plat- 
form/bearing roots over the damping in the autocollimator loops, primarily 
due to the absence of the autocollimator filter. 

The input-output relationships of the gyro stabilized system 
are quite similar to those of the autocollimator stabilized system and 
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Table 5 


AUTOCOLLIMATOR STABILIZATION CONTROL LOOP RESPONSE 
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Fig. 26. MAGNITUDE OF SENSOR RESPONSE TO 




NOISE n -- AUTOCOLLIMATOR STABILIZATION 








Fig. 27. X-Y AXES GSTRO STABILIZATION LOOPS. Only at low frequencies is the viscous cross- 
axis coupling important. 















consequently are not tabulated here. They may be found by simple modi- 
fication of Table 5 , using the new values of the low frequency closed 
loop poles listed in Table 6 . 


Table 6 

GYRO STABILIZATION LOW FREQUENCY CLOSED LOOP ROOTS 


Parameter 

X 

Y 

Units 


5 

5 

rad/sec 





GJ 

C 4 

160 

160 

rad/sec 





0) 

P 

43 

60 

rad/sec 


0.24 

0.27 

— 


Again, only the primary X and Y axes are of later interest. The 
autocollimator filter poles appear only in the autocollimator output 
transfer functions and do not affect the closed loop roots. 



FIG. 28 X-Y AXES GYRO STABILIZATION LOOP ROOT LOCUS 
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Chapter IV 


SPECTRAL ANALYSIS TECHNIQUES 


A. Introduction 

In evaluating the performance of precision instruments, the primary 
interest is in characterizing the errors present in the output signal. 
These errors can be categorized as systematic, random, or input dependent. 

Systematic errors, those which can be described by an explicit model, 
are important to identify since they indicate error sources inherent in 
the design and those sensitive to the environment. Knowledge of the in- 
herent error sources (e.g. , gyro fixed drift) is necessary to verify that 
design specifications are being met, and may lead to performance improve- 
ment through compensation or design changes. It is also important to 
determine the errors due to the laboratory environment (e.g., gravity, 
temperature) to predict performance in the intended application of the 
instrument. 

Random errors can only be described in terms of their statistics. 
These errors cannot be compensated directly; however, the statistical 
knowledge is often used to improve overall system performance through 
filtering or other means of noise rejection. 

In addition to systematic and random errors, which are normally 
attributed to the instrument, a third source of error is uncertainty in 
knowledge of the input to the instrument. For example, case motion of 
a test gyro due to unknown seismic disturbances results in output signals 
which appear to be noise. This noise must be considered random since the 
disturbances cannot be predicted; however, the error could be compensated 
if the disturbance environment were known or measurable. 

The immediate goals of testing are to separate systematic errors 
from random errors, to quantify the errors, and to reduce errors intro- 
duced by the test itself as far as possible. Analysis of the spectral 
content of signals provides considerable insight into noise mechanisms 
as will be shown subsequently. Correlation between signals to separate 
systematic from random errors can also be accomplished in the frequency 
domain. 
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A useful application of spectral analysis is in determining environ- 
mental effects such as vibration, temperature, and gravity. Most instru- 
ments designed to make these measurements are frequency independent only 
within a particular frequency band. For instance, bubble levels, used 
to measure the direction of gravity, indicate with fidelity only at low 
frequencies, whereas seismometers are insensitive at low frequency but 
provide good mid-frequency information. By comparing the measurements 
of such auxiliary devices with the instrument under test and relying on 
each auxiliary device only in its useful frequency band, accurate infor- 
mation about the instrument errors can be obtained across a wide range 
of frequencies. 

The theory behind spectral analysis is well developed [W4], [J2]. 
However, noise is characteristic to each situation and must be analysed 
separately. Only recently has it been practical to use this theory to 
determine specific models, primarily through the development of high 
speed computers and efficient algorithms to compute the spectral functions. 
In Chapter V, the application of spectral analysis to the Fixed Base Simu- 
lator is discussed. To provide a background for this discussion, a brief 
review of spectral theory in a form suitable for digital computation is 
included here. Much of this material in greater detail can also be found 
in [B2] and [K2], The major concepts are first discussed, and then ex- 
amples provided to illustrate the concepts and demonstrate their utility. 

B. Discrete Fourier Transform 

The Fourier transform of a continuous function x(t) is defined as 


X(f) 


■L 


-j2jtft 

x(t) e dt 


and, conversely, the inverse transform is defined as 


x(t) 


-£ 


X(f) e j2 * ft df . 


(4.1) 


(4.2) 
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Digital computation of these transforms require that the data appear as 
discrete points over a finite time interval T. Although uniform time 
spacing of the points is not necessary, the discrete points are assumed 
to be obtained by sampling the continuous signal x(t) every At = 1/f 

s 

sec, where f is the sampling frequency in Hz. The number of data 
s 

points generated is therefore N = T/At. With the data structured in 
this manner, a discrete approximation to the integral in (4. 1) is 


X(f) 


= At 


N-l 


I 

k=0 


x(kAt) e 


-j2jrfkAt 


The overscript C ) is used to denote an approximate or estimated value. 

This function is periodic at the sampling frequency f = 1/At; i.e. 

s 

X(f ± if ) = X(f) i = 0, ±1, ±2 ... (4.3) 

s 

For computational convenience, X(f) is computed only at integer 
multiples of the fundamental frequency, Af = 1/T, of the truncated 
time signal; i.e. 


X(nAf) = At Y x(kAt)e* J23tnAf ^ t . 
k=0 

Since Af = 1/NAt, (4.4) can be rewritten more simply as 


(4.4) 


where 


X(n) 


= At 


N-l 

I 


k=0 


x(k) W nk , 


(4.5) 


W , e- M . 

An approximation to the inverse Fourier transform is 


(4.6) 
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(4.7) 


x(k) = Af y X(n) w' 
n=0 

A 

Both X(n) and x(k) have period N due to the periodicity of W X . 

The discrete approximation (4.5) and (4.7) are related to the Dis- 
crete Fourier Transform (DFT) and the Inverse Discrete Fourier Transform 
( IDFT) treated extensively in [Cl] by 

DFT[x(t) ] = X(n)/Mt 

and 

IDFT [X(n) ] = x(k)/At . 

A more rigorous derivation of the DFT from the Fourier integral can be 
found in [C2]. 

C. Autocorrelation Functions 

The autocorrelation function is useful in determining the dependence 
of a signal on itself displaced in time. It is defined as 

R xx (t,T ) = E [x*(t) x(t+T)] , (4.8) 

where x* (t) is the complex conjugate of x(t). 

If, in addition, the process is wide-sense stationary, its mean is 
constant and the autocorrelation is dependent only on the relative time 
displacement t, 


R xx (T) = :E ^*(t) x(t+T)] . (4.9) 

Furthermore, if x(t) is assumed ergodic, this expression is equivalent 
to 

T 

R nW = 2? / x * (t) X(t+T) dt • < 4 - 10 ) 

Note that R xx (°) = x 2 + where x = E(x), = E(x-x) 2 . 


94 



One discrete approximation is found directly by assuming again uniform 
time spacing of the data over an interval T. 

R x X (m) = iTin / x *(k) x(k+m) in = 0,1, N-l 

k=0 

= 0 m < 0, m > N , 

(4.11) 

where mAt = t, and the time spacing At has been omitted from the 
indices for simplicity. An indirect method of computing the autocorrela- 
tion function will be described after a discussion of power spectral den- 
sity functions. 


D. Power Spectral Density Functions 

The power spectral density function* of a stationary random process 
x(t) is defined as the Fourier transform of the autocorrelation function 


P 

xx 


(f) 


f 


R (t) 

XX 


e"* 2 ** 1 dt 


(4.12) 


Alternatively, 
the spectrally 


the power spectrum can be expressed as a time average of 

** 

decomposed power in the process x(t), 



■ 

m 

2 

p (f) = lim e 

1 

f x(t) e-J 2 * ft dt 


XX T->00 

2T 

~-T 





- 


(4.13) 


The terms "power spectrum" and "PSD" are also used. The use of the 
adjective "power" is not generally descriptive but is included as a 
matter of convention. 

•* 

Although the expected value operator in (4.13) is often omitted in 
engineering texts, it is essential that it be included. This is be- . 
cause (1/2 T) | x(t) e“J 2rtft dt|^ is not a consistent estimator of 

P xx (f); i.e. its variance does not go to zero as time tends to infinity. 
(See [B2] p. 82 and [PI] p. 343, also Appendix A). 


95 



This function is an indication of the frequency distribution of power 
in the signal x(t). If the signal is considered to be the superposi- 
tion of an infinite number of sinusoids of varying amplitude and phase, 
then P(f)Af is the power contributed by the sinusoids in the infinitesi- 
mal frequency band Af centered at frequency f. No phase information 
about x(t) exists in the power spectrum since the autocorrelation func- 
tion is independent of shifts in the time origin. 

By using the Fourier inversion formula, a transform pair between the 
power spectrum and autocorrelation function results, 

R (t) = j P (f) e ^ 2jtfT df . (4,14) 

XX J XX 

•WYl 


Setting t = 0, (4.14) yields the total power in the signal 


R (0) 

XX 




P (f) df = E[x(t) 2 ) = x 2 + a 2 
xx x 


(4.15) 


Similarly, integrating the power spectrum over a finite frequency range 
yields the power in that band. 

A discrete power spectrum can be computed directly by approximating 
(4. 13) using (4. 5) 


/s 

P 

XX 


(n) 


At 

N 


N-l 

I 


k=0 


x(k) 



(4.16) 


= Af | X(n) ] 2 . 


(4.17) 


The power spectrum estimate, divided by Af, is also the least-squares 
estimate of the squared magnitude of the n th complex coefficient in a 
truncated Fourier series representation of x. To show this, assume x 
can be represented by a truncated Fourier series 
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j2itnk/N 


(4.18) 


N-l 

x(k) = > a e 

Z n 

n=0 


It is well known that if a set of measurements x are related to a set 
of parameters a by 


nxl nxp pxl 

x = A a , n > p, rank A = p f (4. 19) 


where the dimensions of the matrices are indicated as overscripts, that 

/N /\ T ^ 

the estimate a which minimizes the mean squared error (x-Aa) (x-Aa) 

satisfies the relationship 


Let 


and 


-1 


a = (A* X A) A* T x . 


(4.20) 


x = 


(X 1 X 2 • • • V 


a = [ ai a 2 


1 T 

V ’ 


A = 


1 1 


1 e 


J2*/N 


1 e 


j2jt (l(-l)/N 


j2jt(N-i)/N 


j2jt (N-l) /N 


We have, by the periodicity of e j 2]Tkn / N > that the columns and rows of 


A are orthogonal, resulting in 


*T 

A A = NI , 

where I is the identity matrix. Therefore from (4. 19) 
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(4.21) 


A 

a 


i 

N 



t 


or 


Finally, 


or 


A 

a 


y 1 x(k> e - J2,tkn/N . 

k=0 





(n)/Af , 


(4.22) 


(4.23) 


(4.24) 


as was to be shown. 

An indirect method of computing a discrete autocorrelation function, 
alluded to in the previous section, is to approximate (4.14) using (4.17) 
and (4.7). 


N-l 

R xx (m) = Af \ | X(n) J W . (4.25) 

n=0 

It is often more efficient to compute autocorrelation functions using 
(4.25) rather than using (4.11) directly, due to time savings afforded 
by the efficient DFT algorithms. However, the results obtained by the 
two methods are not identical. The autocorrelation function computed 
using (4.25) is a circular* correlation related to the direct autocorrela- 
tion function (4. 11) by 


R xx (m) _ — x* (n) x(n+m). Values of the indices falling outside 

the summation interval are shifted into the other end of the interval; 
hence, the term "circular." (See [S3]) 
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(4.26) 


^ c N— in in 

R (m) = R (m) + - R (N-l-m) 
xx N xx N xx 


as can be verified by direct substitution. Only for the case R (m) = 0, 
m > (N-l)/2 are the two discrete autocorrelations equal. However, for 
autocorrelations which go to zero rapidly with increasing values of the 
lag index m, the difference is small. in practice, the overlap in 
(4.26) can be eliminated by appending zeroes to the original data or 
blanking out the last half of the data sequence, (see [B2]). 


E. Cross-Correlation Functions 

The time dependence of two stationary random processes x(t) and 
y(t) can be determined through use of the cross-correlation function 


R (t) = E[x*(t) y (t+T) ] , 

xy 


(4.27) 


or, alternatively, 


R (T) = lira E 

xy t->oo 


By a simple change of variables, 


[A C ‘ 


(t) y(t+t) dt 


b 


(4.28) 


R (T) = R (-t) . 
xy yx 


(4.29) 


If R (x) = x y for all t, where x = E[x(t)] and y =E[y(t)], x(t) 
xy 

and y (t) are uncorrelated. The stronger condition of statistical in- 
dependence of x(t) and y(t) implies uncorrelation. However, the 
converse, that uncorrelation implies independence is not generally true, 
except for Gaussian prpcesses. 

An important application of the cross-correlation function occurs 
in the study of linear systems. Let x(t) be a scalar input to a linear 
system with an impulse response h(t), and let y(t) be the output. 
Then, 


y(t) 


-/ 


x(t-cr) h(q) da 


(4.30) 
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Instead of working with the signals directly, it is often more desirable 
to use only their statistical properties. Multiplying both sides of 
(4.30) by x(t-T) and taking the expected value yields 


R ("0 = f R (t-a) h(cf) do . (4.31) 

xy J xx 

—00 

Given the cross-correlation and autocorrelation functions R (x) and 

xy 

R (t) , the impulse z'esponse of a linear system is determined by (4.31) 

XX 

The solution is easily expressed in terms of the cross-spectral density 
function discussed in the next section. 

A discrete approximation to the cross-correlation function can be 
represented similarly to the discrete autocorrelation function (4. 11) 


N-l-m 

V m) = F5 V x * (k > y< k+m > 

k=i) 


= 0 


m — 0,1, . . . , N-l 


m > N, m < 0 


(4.32) 


Discrete Fourier Transforms can also be used to compute a circular cross- 
correlation function analogous to (4.25). This technique is based on 
the relationship between the cross-correlation function and the cross- 
spectral density function. 


F. Cross-Spectral Density Functions 

The cross-spectral density function is defined as the Fourier trans- 
form of the cross-correlation function 

ym CO 

P (f) = R CO e'j 2 * fT dT , (4.33) 

xy xy 

and, conversely, 


R 

xy 


(f) 



P 

xy 


(f) e 


j2jtfx 


dt . 


(4.34) 
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Note that if x and y are uncorrelated, R (t) = 0, which implies 

xy 

from (4.33) that P (f) = 0. 

xy 

For ergodic processes, the cross-spectrum can also be expressed as 


P , 

xy 


(f) 


lim 

T->oo 


1 

2T 



x* (a) y O') e 


-j2irf (7-a) 


dcr d 7 


(4.35) 


The cross-spectrum is useful in determining the linear relation- 
ships between x(t) and y(t). in particular, for a linear system 
with single input x(t) and output y(t), 


P (f) = H(f ) P vv (f) , (4.36) 

xy xx 

where H(f) is the Fourier transform of the system impulse response. 
From (4.36) it can be seen that 


| P (f) | 

l H(f) l - /V) - (4 ’ 37) 

xx 

and that 

^H(f) = <Cp (f) . (4.38) 

A further relationship between the output power spectrum and the cross- 
spectrum is 


P w (f) = P < f > H * * (4.39) 

yy *y 

Combining (4.36) and (4.39) results in the well known equality 

P yy <f) = | H(f ) 1 2 P xx (f) . (4.40) 

It is evident that computing only the input and output power spectra 
yields only magnitude information about H(f ) . The phase relationship 
between x and y as a function of frequency can only be obtained from 
the cross-spectrum. 
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A discrete approximation to the cross-spectrum can be made by apply- 


ing the DFT to (4.35) 

F> xy (n) = Af xtn) Y(n) n = 0, 1, , N-l . (4.41) 

Directly analogous to the circular autocorrelation function (4,25), 
a circular cross-correlation function can be computed by taking the 


N-l 

✓s c c— ' ^ — nm 

R (m) = Af \ P (n) W m = 0, 1, ..., N-l . (4.42) 

xy / . *y 

n=0 


This circular cross-correlation is related to the lagged product cross- 
correlation (4.32) by 

/s c N-m ^ rn 

R (m) = — rr- R (m) + - R (N-m) m = 0, 1, N-l. (4.43) 

xy N xy N yx 

Again, this formula can be verified by substitution. 

G. Coherence Function 

In the study of linear systems, it is important to have a measure 
of how much of the output is due to the known input. An appropriate 
spectral function to provide this information is the coherence function, 
defined by 


Since 


4 (f) 


p (f) 
xy 


P (f) P (f) 

XX yy 


(4.44) 


V f) t < p xx (f) V f > > 


0 < 7 (f) < 1 . 

— xy — 
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Consider the linear system in Fig. 29. 



Fig. 29. LINEAR SYSTEM WITH MEASURE- 
, MENT NOISE. 

Assume the measurement noise n is uncorrelated with the input x. In 
this case, the spectral functions of interest are 


P 

yy 


(f) 


P (f) + I H(f ) I P = P (f) + P (f) , 

nn 1 ' xx nn ss 


P xy (f) = H(f) P xx (f) . 


(4.45) 

(4.46) 


Forming the coherence function, 


7 2 (f) 

xy 


H(f) 


2 2 
P (f) 

xx 


p (f ) r i H(f) i 2 p (f) + p (f) 

xx 1 1 xx nn 


(4.47) 


or 


2 ... 1 

V f) = ITp 


( f )/p (f) 

nn ss 


(4.48) 


The value of the coherence function is therefore controlled by the ratio 

of the noise power spectrum to the signal power spectrum. For a low sig- 

2 

nal to noise ratio P _/P „ « 1, 7 ( f ) = and the output is said to 


ss nn 


xy 


and the 


be incoherent with the input. For P /P » 1, 7 (f) = 1, 

ss nn xy 

input and output are nearly coherent. Note also that with P ^ 0, 

2 nn 

two values of ] H(f) I obtained from (4.37) and (4.40) are not equal. 


the 


The result using the cross-spectrum, (4.37), is unbiased for no unknown 

2 

noise at the input, whereas (4.40) is always biased except when 7 = 1. 

The former expression generally introduces less bias in the estimate of 
the transfer function. 
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A discrete approximation to the coherence function (4.44) can be 
computed directly from the discrete version of the cross-spectrum and 
the power spectra, (4.41) and (4.17) respectively , 


7 <n) = 

xy 


P (n) 

1 xy 

P (n) P (n) 
xx yy 


(4.49) 


An interesting property of the coherence function estimate can be 
demonstrated for the case when the cross-spectrum and power spectra esti- 
mates are based on just one record of data. Then, 



(n) 


2^* /S /s /X* 

Af X (n)Y (n)X(n)Y (n) 
Af 2 X* (n) X (n) Y* (n) Y (n) 


(4. 50) 


independent of the length of the data record. Therefore, to obtain a 
meaningful estimate of the coherence function, it is necessary to esti- 
mate the cross-spectrum and the power spectral independently. This can be 
accomplished by averaging the estimates from several data records before 
dividing as in (4.49) or by averaging the estimates over nearby frequency 
estimates. Other considerations which must be made in order to obtain 
reliable estimates of the various spectral and correlation functions from 
finite sampled data are treated in the Appendix. 


H. Relationship of Spectral Analysis to Linear Estimation 

The concepts introduced in the previous sections have an important 
application in the area of linear estimation and filter theory. The 
problem is: given measurements, y, of a signal, x, with additive noise, 

n, 

y (t) = x ( t) + n(t) , (4.51) 

determine the linear filter h(t,T) such that the estimate of the signal 
in the form 
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h(t,x) y ( t) dt 


(4.52) 


*o (t) 



T. < t < T f , 


minimizes the mean-square error. In other words, for any other linear 

A 

estimate x, 


El (x - x) j 2 < El (x-x) I 2 , (4.53) 

' o '— 1 1 * 

By the projection theorem, the filter can be shown to satisfy the gen- 
eralized Wiener-Hopf equation 


R (t, a) = 
yx 


-/ 


h(t,T) R yy (T,a) dT . 


(4.54) 


The optimum linear filter can therefore be derived knowing only R 


and R 


yy 


No other information is necessary. 


yx 


Wiener first solved this problem [W4] using spectral techniques by 
•considering only stationary processes and developed a causal (realizable) 
filter. To develop a relationship between the optimal filter and the 
coherence function, consider an optimum non- causal filter by letting the 
limits of integration in (4.52) go to infinity. Assuming stationarity of 
the processes. 


x (t) 


■r 


h (t) y(t-T) dT 


-CO < t < 00 


(4.55) 


and 


R (a) = 
yx 


• / 


h(T) r^(o-t) dT 


-oo < 0 < CO 


(4.56) 


The solution to (4.56) is found by simply taking the Fourier transform 
of both sides 


H(f ) 


P 

yx 

p 

yy 


(f) 

(f) 


(4.57) 
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is found from 



(4.58) 


(4.59) 


(4.60) 


(4.61) 


S = £ Pxx^t 1 - v] df • (4 - 62) 

Equation (4.62) therefore relates the mean-square error to the signal 

power spectrum and to the coherency of the measurement. When the signal 

2 

and measurement are fully coherent - 1) i the mean-square error is 

zero. Conversely, when the noise completely dominates (?y X = 0) » the 
optimum filter and hence the estimate are approximately zero, implying 
from (4.62) that the mean-square error is the total power in the signal 
itself. 

The mean-square error computed by (4.62) is the minimum value attain- 
able from any linear filter, and can be used in determining the amount of 
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data and complexity of filter required in specific applications. Further- 
more, if the signal and noise processes are both Gaussian, this error is 
the minimum for any filter, nonlinear or linear. 

Solution of (4.56) for a causal filter (i.e, h(T) =0, t < 0) using 
Wiener’s techniques is generally practical only when the power spectrum 
of the measurement, p , is rational. In this case a more heuristic 

yy 

solution of the causal filter than that given by Wiener can be found in 
[B6] and [VI]. 

Kalman and Bucy [Kl] advanced the theory of linear estimation by 
finding a solution to the Wiener-Hopf equation (4.54) for finite time in- 
terval, nonstationary processes. The formulation of the filter requires 
certain assumptions which have proven to be valid for a wide variety of 
practical problems. The requirement in Wiener theory for knowledge of 
the multivariable cross-correlation functions between the measurements 
and signals is replaced in the Kalman- Bucy equations by a finite dimen- 
sional model of the signal processes of the form 


x(t) = C(t) 

z(t) 

(4.63) 

z(t) = F(t) 

z(t) + G(t) w(t) 

(4.64) 

y (t) = x(t) 

+ n(t) , 

(4.65) 


where 


E[z(t Q )] = E[w(t)] b E[n(t)] = 0 


w(t) 
n (t) 


[w T (t) n T (T)3 


Q(t) 

T 

LT (t) 


T(t)~ 

R(t)_ 


S(t-T) 


Z (t o ) 


T 

z (t ) 

o 


>J = £ Q , E^z(t o ) w T (t)j = E : 


z(t Q ) n (t) 


]-• 


Instead of determining the optimum filter using spectral techniques, 
the Kalman-Bucy theory expresses the filter as the solution to a set of 
ordinary time differential equations. 
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(4.66) 


z = Fz + K(y-Cz) , z(t Q ) = 0 

K = |zC T + GT^R" 1 (4.67) 

Z = FZ + £F T - KRK T + GQG T , £ (t ) = £ . (4.68) 

o o 

Bryson and Ho have shown the relationship between the two filters 

[BIO]. As a particular case, if F, G, C, R, and Q are constant matrices 

and w(t) and n(t) are stationary processes, the filter gain K may 

reach a steady-state value K . Equation (4.66) may then be written as 

s s 

A A 

z = [F - K C]z + K y . (4.69) 

ss ss 

Laplace transforming (4.69), 

z = [SI - (F-K C)] -1 K y (S) , (4.70) 

s s s s 

where I is the identity matrix. Defining 

®(t) =f _1 [SI - (F-K C))" 1 , (4.71) 

ss 

the solution to (4.69) may be expressed as 


z(t) 



$(t-T) K y(T) dT 
ss 


(4.72) 


Therefore, as t becomes large, the effects of the initial conditions 
diminish so that (4.72) can be written as 


z(t) 



<&(t-T) k y(T) 

b b 


dT 


— oo 

Comparing (4.73) with (4.52) and (4.55), it can be seen that 


(4.73) 
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h(T) = $(T) K 


(4. 74) 


ss 


is the causal Wiener filter. 


I. Illustrati on 

As an example of the utility of spectral techniques in determining 
systematic errors, the situation shown in Fig. 30 was digitally simulated. 



Fig. 30. MULTIPLE SIGNAL MEASUREMENTS WITH 
INDEPENDENT NOISE. 

The two random noise sources n and n were generated independently 

* 1 z 

and are sampled versions of white noise sharply bandlimited at 25.6 Hz. 
Sampling frequency corresponds to 51.2 Hz. The variance of each noise 
source was 51.2 which results in a power spectral density of 1, 0 < f < 
25.6. A small sinusoid at 10 Hz of amplitude 1 was taken to be the sig- 
nal (or systematic error) . Forty data records each 10 seconds long 
(Af =0.1 Hz) were processed and averaged for both measurements y^ and 
y 2 - Figure 31a is a typical time history of the measurements. 

The power spectra of the two measurements are shown in Fig. 31b and 
31c. A spike in the power spectra is apparent at 10 Hz. However, the 
statistical variability of the rest of the spectrum places some uncer- 
tainty in interpreting the spike at 10 Hz as a true component. The cross- 
spectrum between the two measurements Fig. 31d shows the independence of 
the noise, since the level of the noise is attenuated by an average fac- 
tor of 6.3 while the systematic peak at 10 Hz is only slightly attenuated — 
mainly due to statistical variability. 
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Fig. 31a. TYPICAL MEASUREMENT TIME HISTORY. 
The small 10 Hz systematic component is 
not evident. 
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Fig. 31. POWER SPECTRA AND CROSS-SPECTRA OF MULTIPLE SIGNAL MEASUREMENTS. Computing the 
cross-spectrum between a signal and itself displaced in time can uncover the presence of 
a systematic components having relatively long correlation times. 
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Fig. 31e. COHERENCE OF MULTIPLE SIGNAL MEASUREMENTS. The relatively high level 
of coherence at 10 Hz indicates that a common sinusoidal source is present in 
both measurements. The small spike at 3 Hz is due to statistical variability. 


The coherence function Fig. 31e is theoretically equal to 1 at 10 Hz 
and zero elsewhere. However, due to the finite record length, leakage, 
and variability in the estimates of the cross-spectrum and power spectra, 
it only has a value of 0.35 at 10 Hz and near zero elsewhere. Therefore 
some judgment has to be used in interpreting the results. Bendat and 
Piersol [B2] discuss these practical considerations and give measures for 
the expected variation in the spectra and coherence function. Appendix A 
also contains a brief discussion. 

The conclusions that may be drawn from these spectral plots are 

1) the measurements are uncorrelated at all frequencies except at 10 Hz; 

2) a common sinusoid at 10 Hz exists in both measurements; and 3) the 
two measurements are moderately coherent at 10 Hz. Further investigation 
into the statistical properties and computational aspects of the estimate 
must be made before attributing the lack of coherency to variability or 
to other effects introduced by finite record length. 
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Chapter V 


EXPERIMENTAL RESULTS 


A. Introduction 

An application of spectral theory to the modeling of orientation 
sensors is described in this chapter. To demonstrate the application 
realistically, spectra obtained from the gyros, autocollimators, and 
other components of the Fixed Base Simulator are studied. Thus, the 
experimental results presented are primarily intended to demonstrate the 
spectral approach to modeling instrument errors rather than to evaluate 
the performance of the simulator. They present a guide to how spectral 
theory will be useful in the evaluation of the Stanford gyro test of 
relativity. 

It is assumed that the system must be controlled or pointed actively 
to keep the sensor outputs within a representative range. Also, it is 
assumed that the orientation sensors have greater accuracy than available 
laboratory references, and, consequently, must be used to maintain their 
own orientation. As a result, the sensors become active elements of the 
pointing loops, with their output characteristics altered by the system 
dynamics. To properly evaluate the inherent instrument errors, system 
effects must be distinguished. This can be accomplished directly through 
measurements with a more accurate reference, or indirectly by modeling 
the system. These two approaches are discussed in turn. 

A reference sufficiently accurate to test the orientation sensors, 
if available, must necessarily be derived from the sensors themselves, 
since other laboratory references are inadequate. For example, if the 
optical reference is known to have much less error than the gyro at low 
frequencies, then the difference in outputs can be attributed to errors 
in the gyro at these frequencies. Conversely, at high frequencies, gyros 
often have less noise than optical references, implying that differences 
in output are optical errors. This concept of using complementary sen- 
sors, each providing the reference at frequencies where its errors are 
small, could not be exploited in the present research since the autocol- 
limator had superior accuracy throughout the frequency band considered 
(0.1 Hz to 200 Hz). Consequently, the indirect modeling approach was used. 


FNG PAGE BLANK NOT FILING 


115 



The concept of modeling has many interpretations, several of which 
should be distinguished. 

. Writing a set of differential equations from an a priori 
understanding of the system. 

• Deducing a set of equation coefficients from empirical data, 
with or without knowledge of the structure of the equations. 

. Determining a set of equations and coefficients, given an 
output time history and known input. 

• Selecting a set of transfer functions which best match 
observed frequency responses of the system. 

• Making plausibility arguments about system driving noise, 
given the output and knowledge of the system structure. 

• Refining a priori models based on output data and plausible 
models of input noise, and developing descriptors both 
numerical and physical to explain model differences. 

In Chapter III, the first two interpretations were used to gain a basic 
understanding of the simulator. The next two fall under the general 
classification of parameter identification, an area which has recently 
attracted widespread interest. The final two senses of modeling are 
those employed in analyzing the experimental results of this chapter. 


B. Open Loop Sensor Measurements 

The power spectra of signals from the primary autocollimator and 
gyro pickoffs with the sensors operating open loop are first presented 
in Fig. 32 and Fig. 33 as the "true" basic spectra of the instruments. 

The test setups used in obtaining the spectra were sufficiently quiet 
that base motion contributions to the signals are negligible. These 
spectra will be used later as bases of comparison for closed loop results. 


1- Test Setup 

The autocollimator spectra were computed from signals recorded 
while the primary autocollimator was separated from the simulator and 
positioned on a massive marble slab. Relative motion between the auto- 
colliraator and mirror at frequencies of interest was estimated to be 
several orders of magnitude less than the noise equivalent angle of the 
autocollimator. The autocollimator output signals are therefore inter- 
preted as being entirely due to inherent instrument errors. 
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Fig. 32a. OPEN-LOOP POWER SPECTRUM OF THE PRIMARY AUTOCOLLIMATOR X-AXIS CONTROL CHANNEL. 
With the autocollimator resting on a massive marble slab, the effects of base motion are 
negligibly small. 
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Fig. 32b. OPEN-LOOP POWER SPECTRUM OF THE PRIMARY AUTOCOLLIMATOR X-AXIS INSTRUMENTATION CHANNEL. 
The primary difference between this spectrum and the previous one is the filtering of noise 
above 40 Hz in the control channel. 
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Fig. 33a. POWER SPECTRUM OF THE X-AXIS INSTRUMENTATION CHANNEL OF X-Y GYRO WITH PLATFORM 
CLAMPED. The concentration of power at discrete frequencies is characteristic of ball 
bearing gyros. 
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Fig. 33b. POWER SPECTRUM OF THE Y-AXIS INSTRUMENTATION CHANNEL OF X-Y GYRO WITH PLATFORM 
CLAMPED. Several of the spikes are also present in the X-axis, indicating a common source. 


The gyro power spectra were obtained with the gyro mounted on 

the platform, but with the stabilization loops disabled and the platform 

rigidly clamped to the rotary table. With this rigid arrangement, any 

motion sensed by the gyro is due primarily to microseisms. Bradner and 

Reichle [B7] state that for an average test site the acceleration power 

2 

spectrum in the range 0.5 to 50 Hz is approximately 0.001 g /Hz. Angular 
motion can be deduced from the linear acceleration knowing the soil shear 
propagation velocity and frequency of the underlying waves. A relation 
given in [FI] is 


a 


5y 

Vf 


where 

O' = peak angle (prad) 
y = peak vertical acceleration (pg) 

V = soil shear propagation velocity (ft/sec) 
f = wave frequency (Hz) 

Using a value of V = 400 ft/sec, the angular power spectral density at 

-7 2 

0.5 Hz is approximately 6 x 10 prad /Hz and varies inversely with the 
square of frequency. This value is considerably below that observed in 
the gyro spectra. Gyro float motion due to mass unbalance and seismic 
accelerations are also negligibly small at frequencies above 0. 1 Hz. 
Vibration induced by the gyro rotor bearings and electronic noise are the 
most probable sources of the pickoff noise. 

2. Graphical Interpretation 

Several general characteristics of these spectra should be 
noted. Firstly, the noisiness or jitter of the plots is due to insuffi- 
cient averaging; only five estimates of each frequency point were averaged 
to obtain the values presented. Later spectral plots will appear smoother 
since the number of estimates in each average was increased. Secondly, 
frequency resolution of these spectra, as well as all other plots presen- 
ted is approximately equal to the smallest frequency plotted. In the 
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present case, resolution is about 0. 1 Hz. Thirdly, low-pass anti-aliasing 
filters with cutoff frequencies between 100 and 150 Hz were used to pre- 
vent spectral distortion below 100 Hz. Dramatic decreases in power above 
the cutoff frequency are evident in all the spectral plots. Finally, the 
spectral plots shown are two-sided, requiring the area under the curve to 
be doubled in computing the actual power in any particular frequency band. 

3. Primary Autocollimator Spectra 

The basic autocollimator spectra of Fig. 32, taken from the X- 
axis control and instrumentation channels, can be separated into three 
frequency bands. In the mid-f requency band, i.e. between 3 Hz and 35 Hz 
for the control channel and between 3 Hz and 80 Hz for the instrumenta- 
tion channel, the spectra have constant spectral density, a characteris- 
tic of white noise. The level of noise in this band was predicted during 
design based on a shot noise model. The noise is almost entirely due to 
dark current in the autocollimator photodetector. Amplifier noise is 
negligibly small. 

Below 3 Hz the power spectra vary inversely with frequency. 

Noise exhibiting this type of spectral behavior is often called 1/f or 
flicker noise [B3]. It can be generated electrically due to surface ef- 
fects of solid state devices. However, since the autocollimator is an 
ac device, electronic 1/f noise is shifted to the 2.2 kHz switching 
frequency. Demodulation and filtering effectively eliminate this noise 
from the output. The 1/f characteristic in the spectra is therefore 
not likely caused by the electronics. A more plausible source is dis- 
tortion of the infra-red rays due to turbulence in the air between the 
autocolliraator and the reference mirror. This explanation was not well 
established. 

Above the mid-frequency band, the noise decreases due to fil- 
tering. The spectral curves in these high frequency bands correspond to 
the squared magnitude of the filter transfer functions multiplied by the 
constant spectral density of the white shot noise source. The response 
is that predicted from a deterministic analysis of the instrument. This 
observation is important since it demonstrates the validity of the statis- 
tical spectral approach. Here the coloration of the noise clearly suggests 
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a model of the instrument without having to inject external signals to 
determine its response. 

The control channel break frequency of 40 Hz corresponds to 
the two-pole output filter in the autocolliraator electronics, while the 
instrumentation channel attenuation is due primarily to the anti-alising 
filter. Flattening of the spectra above 100 Hz is caused by the equiv- 
alent noise generated by quantizing the analog signal in the analog-to- 
digital converter [W3]. The autocollimator rms noise equivalent angle is 
0.55 |_irad, found by taking the square root of twice the area under the 
instrumentation channel spectrum. Similar measurements were used in ob- 
taining the autocollimator performance data presented in Table 1 in Chap- 
ter II. 

4. Gyro Spectra 

The spectra from both axes of the X-Y gyro, shown in Fig. 33, 
are characterized by a white background with numerous spikes at discrete 
frequencies. This type of spectrum is typical of ball bearing gyros. 

The discrete components at 23 and 28 Hz appearing in both X and Y axes 
are not consistently present in spectra taken at other times and are of 

unknown origin. Three other large spikes appear in both X and Y axes 

at 71, 102 and 187 Hz. The 187 Hz discrete component is due to vibration 
at the rotor spin frequency. The induction spin motor field rotates at 

200 Hz and a 13 Hz rotor slip frequency produces the torque required to 

overcome frictional losses. A small spike in the spectra near 13 Hz pos- 
sibly corresponds to this slip frequency. The rotation of the rotor 
bearing retainers is a potential source of a large discrete frequency 
component. Unfortunately, the configuration of the rotor bearings is not 
known so the rotation frequency of the retainers cannot be computed. It 
is known, however, from available gyro documentation that the rotor is 
actually composed of two flywheels. It is likely that the bearing re- 
tainers of these flywheels do not have exactly the same rotation frequency 
due to manufacturing tolerances. The apparent large spike at 71 Hz is 
actually composed of two large spikes separated by about 0.8 Hz. Beating 
in the gyro pickoff signal, which would result from two nearby frequency 
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components, was observed in strip chart recordings of early simulator 
runs, although this phenomenon was not apparent in later tests. The com- 
puted spectra, however, consistently show two distinct frequency compo- 
nents, strongly indicating that rotation of the bearing retainers is re- 
sponsible for the concentration of signal power near 71 Hz. Another 
possible, though less likely, source is nutation of the float, which 
occurs at a natural frequency h/l ^ = 625 rad/sec. Inaccuracies in the 
measurement of the float inertia If. could account for the discrepancy 
in frequency. The low damping ratio (£ = 0.07) of the nutation mode 
would allow large increases in power near its natural frequency. 

The very large spike at 102 Hz is of unknown origin. Nutation 
is a likely source because of the close agreement in frequency. However, 
the height of the spike implies that a nutation damping ratio less than 
0.02 is required, which is substantially less than the 0.07 value esti- 
mated. Other possible sources are electronic noise, irregularities in 
the rotor bearings, or a mechanical resonance. 

It should be observed at this point that the power spectra only 
suggest noise mechanisms, and that further experiments are usually neces- 
sary to verify any particular model. The explanation given above for the 
discrete concentrations of signal power need further substantiation; they 
are, however, plausible based on the available evidence. 

5. Azimuth Autocollimator Spectra 

The power spectrum of the azimuth (Z) reference autocollima- 
tor, Pig, 34, is a clear example of 1/f noise. The data were obtained 
from an experimental setup similar to that used for the primary autocol- 
limator. The azimuth autocollimator is a dc device, and consequently is 
more sensitive to low frequency noise than the primary autocollimator. 

As evident in the spectrum, the output signal is composed entirely of 
1/f noise, at least out to the anti-aliasing filter break frequency. 

Modeling of 1/f noise is especially difficult since it can- 
not be described using a finite number of lumped parameters. In other 
words, a state variable model would require an infinite number of equa- 
tions to specify the noise. 
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An approximate time domain solution to this problem is given 
in [Tl], The noise is treated as the sum of outputs of first order lag 
networks driven by independent white noise sources. The lag time con- 
stants must be distributed throughout the frequency band of interest. 
Accuracy of this approximate model is determined by the number and dis- 
tribution of networks. 


C. 


Coherent Measurements 


It is important to be able to distinguish platform motion from the 
inherent errors of the orientation sensors. With ideal instruments, 
corresponding channels should have the same output: that due to platform 
motion. The degree to which signals are different therefore depends on 
the amount of noise or instrument error present. As discussed in Chapter 
IV, an excellent measure of the amount of noise present in instrument 
outputs is the coherency function. 

An experiment was conducted to demonstrate the use of this function. 
The intention was to raise the signal well above the noise in a particu- 
lar frequency band. To do this, the platform gimbal torquers were excited 
by a low frequency white noise source band limited to 30 Hz while the 
simulator was autocollimator stabilized. Referring to Fig. 24, the ex- 
citation signal was introduced into the X axis as V T . Its power spec- 
trum is shown in Fig. 35; the total power corresponds to an rms excita- 
tion torque of 0. 0078 N-m, which is "small" but larger than the noise. DC 
excitation and response were removed from the computation of the spectrum 
in the manner discussed in Appendix A. 

The motion induced by the excitation torque was sensed by the gyro, 
and the output spectrum computed, Fig. 36. Angular rms amplitude is ap- 
proximately 13 prad. The coherence between the excitation voltage and 
the gyro signal is shown -also in Fig. 37. Coherence is excellent up to 
and slightly beyond the stabilization loop bandwidth, approaching unity 
at these low frequencies. The gyro output below the servo bandwidth is 
therefore almost entirely due to the excitation. At very low frequencies, 
the coherence is not as high, indicating that gyro drift torques or other 
disturbances are influencing the gyro output. Above 20 Hz, the coherence 
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Fig. 35. POWER SPECTRUM OF EXCITATION SIGNAL V T INTRODUCED INTO X-AXIS AUTOCOLLIMATOR STABILIZATION 
LOOP. Attenuation of the signal above 30 Hz is due to filtering in the low-frequency noise generator. 
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Fig. 36. POWER SPECTRUM OF THE GYRO X-AXIS DUE TO EXCITATION V Td> Power below 15 Hz is due to 
base motion while high frequency power is instrument error. 



is nearly zero, indicating that other gyro noise or disturbances are 
responsible for the output. This high frequency portion" of the gyro 
spectrum is nearly identical to the open-loop gyro spectrum. Fig. 33a. 

The similarity implies that pickoff noise is the source of the high fre- 
quency power in both the open-loop and torquer- excited spectra. 

In addition to determining the coherence between the input excita- 
tion and the output, a comparison of the outputs of the gyro and auto- 
collimator was also made. Figure 38 is the power spectrum of the auto- 
collimator. It is quite similar to the gyro spectrum except at high 
frequencies, a fact further substantiated by the coherence between the 
gyro and autocollimator, Fig. 39. The high coherence at low frequencies 
indicates that platform motion is primarily responsible for the output 
of the instruments. Instrument errors at these frequencies are negligible 
by comparison. The reduction in coherence below 1.5 Hz is similar to that 
seen in the previous coherence function, Fig. 37. At higher frequencies, 
gyro pickoff noise causes the two signals to be incoherent. The sharp 
reduction in coherence above the servo bandwidth is evidence of the high 
level of discrimination afforded by the coherence function. 


D. Transfer Function Modeling 

It is not always possible or desirable to excite a system in which 
precision orientation sensors are being evaluated. The sensors may be 
driven out of range, or the excitation magnitude may obscure small but 
important error sources. Furthermore, the excitation may actually alter 
the physical situation so that performance measurements are no longer 
valid. Consequently, it is desirable to model system transfer functions 
using only the noise inherent in the system. The transfer functions may 
then be used in determining the effects of the system on the sensors. 

Figure 40 is the amplitude ratio of the transfer function of the X- 
axis compensation network computed from measurements of the network input 
and output. The simulator was autocollimator stabilized. No external 
excitation was used; the input to the network being the output of auto- 
collimator control channel plus a small but unknown amount of electrical 
pickup noise. The transfer function derived in Chapter ill by a network 
analysis, Fig. 22, is superposed for comparison. 


These two methods of 
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Fig. 38. POWER SPECTRUM OF THE AUTOCOLLIMATOR X CONTROL CHANNEL DUE TO EXCITATION V T . Low- 
frequency attenuation is due to the integral control feedback of the compensation network, 
while high frequency attenuation results from platform inertia and electronic filtering. 
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Fig. 40. COMPARISON OF EXPERIMENTAL AND DESIGN COMPENSATION NETWORK TRANSFER FUNCTIONS 
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Fig. 41. COHERENCE OF COMPENSATION NETWORK INPUT AND OUTPUT. A probable source of the coherent 
components at 200 Hz is electrical pickup from the gyro spin motor drive electronics. 



obtaining the transfer function agree to within a factor of two. This 
is considered good agreement for the amount of data processed and the 
amount of internal noise picked up from surrounding circuits. The mod- 
erate level of coherence in the coherence function, Fig. 41, indicates 
the presence of a significant amount of pickup noise. 

This example of computing transfer functions was selected since the 
actual transfer function was accurately known and could be used as a 
basis of comparison. The technique, however, is applicable in determining 
input-output relationships throughout the system. 

E. Closed Loop Sensor Measurements 

The power spectra of the sensors, mounted on the simulator and op- 
erating as integral parts of the stabilization loops were calculated. 
Comparing these spectra, shown in Figs. 42 and 43, with the open loop 
spectra of Figs. 32 and 33 clearly demonstrates the need for an accurate 
system model in order to deduce sensor performance from closed loop data. 

1. Autocollimator 

The power spectrum of the autocollimator X control channel, 

Fig. 42, was computed from data taken with the simulator in Mode 3 (auto- 
collimator stabilization). This spectrum is to be compared with the open 
loop spectrum. Fig. 32a.' Only above 40 Hz, where the attenuation of the 
output filter is evident, are the spectra comparable in form. In this 
band the level of the closed loop plot is ten times greater than that of 
the open loop. The increased level of noise over the open loop measure- 
ment is most likely due to electrical pickup noise in the much harsher 
simulator electrical environment. Below 40 Hz almost no similarity exists 
because the operating conditions are completely different. Differences 
in the two plots must be explained by considering overall system behavior. 
In order to do this, the simulator model developed in Chapter III is used. 
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Fig. 42. POWER SPECTRUM OF THE X-AXIS CONTROL CHANNEL OF THE PRIMARY AUTOCOLLIMATOR IN MODE 3 
(AUTOCOLLIMATOR STABILIZATION) . 


Referring to the autocollimator stabilization loop block diagram, 
Fig. 24, it can be seen that numerous paths exist between the various dis- 
turbance inputs and the sensor output. Assuming the simulator is ade- 
quately described by the model, each output is a combined response to all 
the disturbances. The responses have been evaluated previously in Table 
5 in terms of the closed loop transfer functions. Since the system was 

not commanded (7 = V_ = 0) , the most likely causes of the autocollima- 

c T d 

tor spectrum are the autocollimator error n , disturbance torques M , 

A 

or amplifier noise, also considered as a disturbance torque. Mechanical 
vibration is another possible source. 

As a starting point in developing a model which explains the 
autocollimator spectrum, the pronounced decrease in power at 3.5 Hz is 
first noted. An explanation of this dip follows. 

Within the bandwidth of the system (15 Hz) the closed loop trans- 
fer function between an input and the output is approximately 1/H(s) , 
where H(s) is the transfer function of the feedback path. A dip in the 
transfer function within the system bandwidth must therefore be caused by 
a peak in H(s). Although resonant peaks can arise from several sources 
(e.g. gyro, reference mirror mount, autocollimator mount, flexible modes), 
the most likely source is a platform gimbal bearing resonance. Other res- 
onances tend to be much higher in frequency than the observed 3.5 Hz, 
while a soft elastic restraint in the platform gimbal bearings coupled 
with the large gimbal inertia could produce the low frequency resonance. 
Referring to Fig. 24, the only noise input which has this potential source 
in the feedback path is the autocollimator noise n . 

A 

Based on this observation, an investigation of bearing behavior 
was undertaken. The study is described in Section F of this chapter. The 
main result, already accounted for in the simulator model of Chapter III, 
is that the bearings act like springs for small motion. 

Using this model, the closed loop transfer function between n 

A 

and the autocollimator control channel output, 7 , is expressed in 

Table 5 of Chapter III as ( • The observed spectrum is sim- 
ilar in form to the magnitude plot of this transfer function shown in Fig. 
26a; the only differences occuring at low frequencies where the observed 
spectrum is higher and decreases only slightly with decreasing frequency. 
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Therefore, with a high level of confidence, the decrease in power at 
3.5 Hz can be attributed to a platform gimbal bearing resonance, and the 
entire autocollimator spectrum attributed to the response of the stabili- 
zation loop to the autocollimator noise spectrum of Fig. 32a. 


2. Gyro 

The spectrum computed from the simultaneously recorded X-axis 
gyro pickoff signal, Fig. 43, is also consistent with these results. The 
model response of the gyro to white autocollimator noise has been shown 
previously in Fig. 26c. Above 15 Hz the gyro spectrum is essentially the 
same as the open-loop results obtained with the platform clamped. Below 
15 Hz the closed-loop spectrum is a combination of pickoff noise and plat- 
form response to autocollimator noise. In the band around 9 Hz, coherence 
exists between the gyro and autocollimator signals. This can be seen in 
Fig. 44. Platform motion is significantly larger than the autocollimator 
noise equivalent angle near 9 Hz due to peaking of the lightly damped 
stabilization loop. The signal to noise ratio in both sensors is there- 
fore increased, a condition reflected by the relatively high level of 
coherence. The signals are also moderately coherent at low frequencies 
for unknown reasons. 

The gyro spectrum can also be used to deduce the level of motion 
of the platform. Again from Fig. 26c, the gyro (and platform) do not re- 
spond significantly to autocollimator noise above 15 Hz. Since the gyro 
spectrum below 15 Hz is primarily due to platform motion induced by auto- 
collimator noise, the total power in this band is a reasonable estimate 
of the variance of the platform angle. This excludes any motion of the 
reference mirror at frequencies less than 0.4 Hz, which should be small 
due to the rigidity of the autocollimator reference mirror mount. The 
computed rms platform angle at frequencies greater than 0.4 Hz is approxi- 
mately 1.6 jirad. Although no data were obtained below 0.4 Hz, very low 
frequency motion due to autocollimator noise should not be significant 
because of the excellent null stability of the autocollimator. Mechani- 
cal creep and thermally induced distortions are more likely sources of 
low frequency error. 
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Fig. 44. COHERENCE OF THE X-AXIS AUTOCOLLIMATOR CONTROL AND GYRO INSTRUMENTATION CHANNELS. 



Corresponding Mode 3 plots for the Y axis are shown in Figs. 
45, 46, and 47. The behavior is similar to the X axis, the most notable 
differences being the general shift toward higher frequencies of the main 
features, and the broader band of the dip in the autocollimator spectrum. 
The large spike at 200 Hz in the autocollimator spectrum is probably due 
to interference at the gyro spin motor drive frequency. The results 
are otherwise consistent with the simulator model developed independently 
in Chapter III. The computed rms platform angle is 3.6 |jrad. 

3. Cross- Axis Coherence 

The coherence between the X and Y autocollimator control 
channels was computed and is shown in Fig. 48. The high level of coher- 
ence above the stabilization loop bandwidth of 15 Hz indicates that the 
two channels have a common noise source. As discussed previously, this 
source is the dark current in the photodetector. At low frequencies the 
stabilization loops reduce the response to the noise to a sufficiently 
low level that other noise sources become important. Consequently, the 
coherence is also small. 

The analogous coherence function of the two X-Y gyro instru- 
mentation channels is shown in Fig. 49. The two channels are nearly in- 
coherent everywhere except at a few discrete frequencies. This indicates 
that most of the noise arises in the separate pickoff electronics rather 
than from a common signal generator excitation noise or a platform motion 
about an axis skewed with respect to both X and Y. 


F. Platform Gimbal Bearing Restraint 

The decrease in power in the 3. 5 - 10 Hz band consistently observed 
in the closed-loop autocollimator spectra was not anticipated. An orig- 
inal model for the platform gimbal dynamics consisted of the gimbal iner- 
tia and a small amount of bearing friction. For the small platform motion 
present, it was not known whether a viscous or coulomb (velocity indepen- 
dent) friction model was appropriate. Using a viscous model, the observed 
decrease in power could not be predicted. On the other hand, it was not 
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Fig. 46. POWER SPECTRUM OF THE Y-AXIS INSTRUMENTATION CHANNEL OF THE X-Y GYRO IN MODE 3. 


1*0 


10*0 

Frequency <Hz) 


100*0 


Fig. 47. COHERENCE OF THE Y-AXIS AUTOCOLLIMATOR CONTROL AND GYRO INSTRUMENTATION CHANNELS 
both e instr^ents? tl0n ^ ^ ^ Stabilization loo P cross-over frequency is sensed * 



Coherence 


l.o 


0*8 


0 • 6 J 


0*4 J 


0*2 


0*0 J 


l r i > i — r — i — i — r-» 

0*1 1*0 


1 1 — • — i i i r i i i — i — i — i — » t i 

10*0 100*0 


1 — I — 1 — I I II 
1000*0 


Frequency (Hz) 


Fig. 48. COHERENCE OF X AND Y AUTOCOLLIMATOR CONTROL CHANNELS IN MODE 3. 
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Fig. 49. COHERENCE OF X AND Y GYRO INSTRUMENTATION CHANNELS OPERATING IN MODE 3. 


known how coulomb friction would manifest itself in the sensor power spec- 
trum. A further indication that the initial platform model was not ade- 
quate was the increase in frequencies of the observed spectral peaks over 
those predicted from the linear model. 


1 • Identification and Measureme nt 

In order to isolate the mechanism by which power is apparently 
absorbed, a study of the transfer functions of the individual system com- 
ponents was made. Transfer functions were determined using the spectral 
techniques already described, and also by inserting sinusoidal voltages 
from a function generator into the autocollimator and power amplifier 
command inputs. All computed transfer functions except the platform dy- 
namics agreed well with the previously established model. 

The platform response was determined by commanding the gimbals 
to oscillate sinusoidally at a very low frequency and simultaneously 
recording the autocollimator instrumentation output, the power amplifier 
input, and the voltage across the gimbal motor armature. The gimbal motor 
armature voltages are proportional to the motor current, and thereby the 
torque, since back EMF of the gimbal motor is negligibly small. The com- 
mand was summed into the autocollimator control channel output with the 
simulator in Mode 3. The instrumentation channel indicated the actual 
gimbal angle. Figure 50 is a representative low frequency time history 
of the recorded signals. It can be seen that as the gimbal angle reaches 
its maximum displacement and starts returning to null, the applied torque 
unloads rapidly and then steadily increases until the velocity of the 
gimbal again changes polarity. 

A cross plot of platform angle and torque taken at low frequency 
(0.4 Hz), Fig. 51, illustrates the hysteresis characteristic of coulomb 
friction. Identical behavior was observed at lower frequencies. The power 
amplifier input was used as the measurement of applied torque since it is 
quieter than the motor voltage. A large amplitude oscillation, 1000 prad 
peak to peak, was used to increase the signal-to-noise ratio. The ampli- 
tude of motion (500 prad) is more than two orders of magnitude greater 
than that experienced during normal simulator operation, and is not repre- 
sentative of small motion behavior. 
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Fig. 50. REPRESENTATIVE LOW FREQUENCY TORQUE— ANGLE TIME 
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Fig. 51. LARGE AMPLITUDE, LOW- FREQUENCY BEARING 
TORQUE-DISPLACEMENT RELATIONSHIP. 



Fig. 52. FREQUENCY DEPENDENCE OF RESPONSE OF 
PLATFORM GIMBAL ANGLE TO MOTOR TORQUE. 
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The ratio of large amplitude motion to the peak applied torque 
taken at different frequencies, plotted in Fig. 52, is typical of a spring- 
mass system. However, the constant ratio at low frequencies can also be 
explained by coulomb friction. The ratio is amplitude dependent and not 
frequency dependent. At high frequencies where the inertia of the gim- 
bal dominates, the ratio decreases as the square of frequency. Signals 
recorded in the transition region between 1 and 3 Hz were erratic in form, 
with time histories a combination of those observed at both low and high 
frequencies. 

Two other low frequency torque-angle cross plots were obtained: 
one at moderate amplitude, the other at very low amplitude, Fig. 53 and 
Fig. 54. In both plots the overall behavior has changed from primarily 
hysteresis to one better described as a spring. Hysteresis is still 
present, but to a much less extent. Data obtained at very low amplitude 
were heavily corrupted by noise and jitter, and some visual averaging was 
necessary to read the values plotted. Results in this case were consis- 
tent, however, as were those at all amplitude levels. 

2. Related Experiments 

The behavior just described can be attributed to the platform 
gimbalbearings . Similar results have been obtained in tests conducted 
by Aerospace Corporation [Dl] and Ball Brother Research Corporation (BBRC) 
[C7] under different conditions and with larger ball bearings. Aerospace 
considered three spherical bearings on a flat surface loaded by another 
plate. Dynamical response to small displacements was observed. BBRC 
developed a fixture in which the outer bearing race is attached to a worn 
wheel driven by a geared dc motor. Displacement of the inner race is 
measured by a linear variable differential transformer (LVDT) . Variable 
bearing preloads can be established by suspending weights from the inner 
race. The inner race is driven by a precision dc motor. With the outer 
race fixed, calibrated torques can be applied by the precision motor and 
the inner race displacement recorded to obtain static friction measure- 
ments. Alternatively, running friction is measured with the outer race 
driven at a known rate and the LVDT output fed back to the motor through 
servo electronics to hold the inner race at null. The motor torque 
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Fig. 53. MEDIUM AMPLITUDE BEARING TORQUE- 
DISPLACEMENT RELATIONSHIP. 



Fig. 54. SMALL AMPLITUDE BEARING TORQUE- 
DISPLACEMENT RELATIONSHIP. 
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required is then just equal to the friction. 

Primary results of the BBRC tests were: 

1. Displacement is proportional to torque at low torque levels. 

2. The bearings do not completely return to their original 
position after torque is applied and then removed. This 
is indicative of losses due to coulomb friction. 

3. Displacement increases indefinitely with time above a 
certain minimum torque level. 

4. The minimum breakaway torque level increases with increasing 
preload. 

3. Physical Interpretation 

Dahl [Dl] has developed a solid (or coulomb) friction model which 
has been used in a simulation to predict bearing behavior accurately. A 
summary of Dahl's physical explanation of bearing behavior is given below. 

Bonds are formed between the bearings and their races, with the 
strength of the bond increasing with preload. As torque is applied, the 
bearings tend to rotate, placing the bonds in tension in the rear contact 
area and in compression in the leading contact area. Figure 55 illustrates 
this concept. 



Fig. 55. BEARING STRESSES UNDER APPLIED TORQUE. 
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For sufficiently low torque levels, the contact stresses are 
proportional to displacement and to the applied torque. As torque is 
increased, permanent deformation of the interface first occurs, followed 
at higher levels by rupturing of the trailing edge of the contact area. 
Depending on geometry and surface properties, at some point the maximum 
restoring torque of the bearing is reached. Further increases in torque 
tend to accelerate the bearing. 

4. Mathematical Model 

In developing a mathematical model of rolling friction, Dahl 
first expressed the friction as the solution of the differential equation 


dT dT • 
dt = d0 y 


(5.1) 


Therefore, with known initial conditions, and a model of the relationship 
between torque and angular displacement, dT/ d 0, the frictional torque 
can be specified. Based on data obtained in the Aerospace study discussed 
above, a model was postulated, 


dT 

d0 

= n (T-T ) 2 
* o 

(0 > 0) 

(5.2a) 

dT 

, v 2 



d0 

= i<t+V 

(0 < 0) 

(5.2b) 


where r\ is a constant, and T q is the maximum frictional torque that 
can be produced (running friction). Excellent agreement between computed 
and experimental values of friction was obtained by employing this square- 
law model in a computer simulation. 

To compare the model with the results obtained from the Fixed 
Base Simulator, (5.2a) and (5.2b) are solved: 


0 + c, = 


riT (1-T/T ) 
1 o o 


(0 > 0 ) 


(5.3a) 
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-1 


(0 < o) 


(5. 3b) 


6 + c 


2 


nT (1+T/T ) 
1 o o 


where c and c are constants. 

J. ^ 

For small motion, T « T , which starts at 6 = 0 and T = 0, 

o 

both Eqs. (5.3a) and (5.3b) can be expanded to yield, approximately, 



This relationship is characteristic of a spring with elastic coefficient, 
2 1 

•nV For large motion, 6 » — — — , the magnitude of friction approaches 

o 

the constant value T , characteristic of coulomb friction. 

o 

Equations (5.3), normalized as discussed subsequently, are 
plotted in Fig. 56. The angle is varied sinusoidally so that the model 
plots can be compared with the experimental frictional torque data of 
Figs. 51, 53, and 54. Several amplitudes of motion are displayed to show 
the strong amplitude dependence of the frictional torque. The location 
of the origin is not of interest in displaying the amplitude dependence, 
so the constants of integration, c and c , were selected to keep the 
figure centered. 

The two constants and the maximum frictional torque, T , de- 
ni 

veloped for any particular amplitude, G , can be found by simultaneous 
solution of the three equations 



+ c. 



6 

o 


+ C. 


h T , 


1 - 




(5.4a) 


(5.4b) 


(5.4c) 
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Fig. 56. SQUARE-LAW MODEL RELATING BEARING FRICTION 
TORQUE TO ROTATIONAL DISPLACEMENT. 
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The solutions are 


T 

m 

T~ 
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2t]T 0 
1 o o 


Jl + 


<2 r\ t e y 

'o o o 


1 


(5. 5a) 



(5.5b) 


(5.5c) 


It is convenient to normalize the angle with respect to the 

non-dimens inal parameter HT : i.e. 

o 


6 = t)T q 0 (5.6a) 

®o =1 » T o S o’ (5 ‘ 6b) 


and to normalize the frictional torque with respect to the maximum pos- 
sible or running torque T • 

o’ 


T 



(5.7a) 



(5.7b) 


This particular normalization does not affect the shape of Figure 56. 
The normalized solutions to Eqs. (5.3a) and (5.3b) then become 


T 


T 2 + (1-T 2 )0 
m m 

1 + (1-T 2 )0 
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0 > 0 
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_o _2 

T - (1-T )0 
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# -2 - 
i - (i-T ye 
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(5. 8a) 


(5.8b) 
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where 


m 


■ ■ *] 


(5.9) 


The basic parameters determined by fitting the model to the 
experimental data of Figs. 51 , 53 , and 54 , are 


T = 0.016 N-m 
o 

0 

T) = 1.6 X 10 (N-m) 


5. Describing Function 

A s'inusoidal input describing function [G2] for the nonlinear 
torque has been derived. from Dahl's model. Using the symmetry indicated 
in Fig. 56, the in-phase and out-of-phase components of the describing 
functions are calculated from 


—2 2 — 

r 2it T + (1-T )9 cos \|r 

T cos 0 = - I — „ m _ ° cos \Jr d\Jr (5. 10) 

* A 1 + (1-T )6 cos \Jr 
m o 


and 

2 r 2lt T 2 + (i-T 2 ) 0 sin ^ 

T sin Cp = — / — _ 2 m _ ° sin di|r (5. 11) 

* “at 1 + (1-T )0 sin \|r 
m o 


where = T D / T 0 is the amplitude, and cp is the phase angle of the 
describing function. Straightforward evaluation of (5.10) and (5.11) 
leads to 


T 

6 
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e 

o 




9 , i + d-T 2 )e 

2 £ n f m ° 


ate \1 - (1-T 2 ) ~ 
o \ m o 


(5.12) 
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This describing function is displayed in Fig. 57. 

PHASE, £ 



Fig. 57. SINUSOIDAL INPUT DESCRIBING FUNCTION OF A SQUARE- 
LAW BEARING FRICTION MODEL. 


For 6 ^ « 1, (5.12) can be simplified to 


_D 0 j0 ~ ^ (spring) 

Q 

o 

Likewise, for 0 q » 1, 


— ^ = j (coulomb) 

e *6 

o o 
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Chapter VI 
CONCLUSION 


A. Summary 

Analyzing the signals from precise orientation sensors using modern 
digital spectral techniques is an effective means of identifying error 
sources and modeling their effects. However, there are numerous consid- 
erations in converting raw spectral information into meaningful inter- 
pretations of instrument errors. The dynamical behavior of the system 
in which the sensors are operating can have a profound effect on the 
output spectra of the instruments. To interpret experimental results 
properly, it is essential to have an accurate model of the system, in 
terms of its natural behavior and disturbance environment. Operating the 
sensors as part of a closed stabilization loop requires that the paths 
through which instrument errors propagate before being measured be well 
understood . 

Since instrument errors are created by both internal sources and 
by system disturbances, it is necessary to separate the sources to identify 
the inherent instrument errors. The coherence function was demonstrated 
to be useful in accomplishing this task. A further useful approach is 
to make plausible arguments about the disturbance spectra and compare 
the experimental results with those obtained by analysis of a system 
model. This technique was successful in discovering the presence of the 
unanticipated elastic restraint properties of the simulator platform 
gimbal bearings. 

Once discovered, further direct measurements of the bearing be- 
havior were made. Elastic restraint coefficients of approximately 400 
N-m/rad were observed, with coulomb friction becoming dominant at 
rotation displacements greater than 100 prad. The relationship between 
frictional torque and angular motion agreed well with a square-law 
model of rolling friction discussed in Chapter V. The model predicts 
spring-like behavior with low amplitude motion and coulomb friction at 
large amplitudes. A sinusoidal-input describing function of the fric- 
tional torque, derived from the model, indicates the strong amplitude 
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dependence and frequency independence of bearing friction. 

The pointing accuracy of the Fixed Base Simulator, developed as an 
analog to the proposed laboratory test of the Stanford Relativity Experi- 
ment satellite, is limited primarily by the noise in the orientation 
sensors. Other disturbances are negligible by comparison. The observed 
spectra of the sensors operating as integral parts of the stabilization 
loops can be attributed to sensor noise as modified by the loop dynamics. 
Platform pointing errors in the band 0.4 to 200 Hz with the primary 
autocollimator serving as the controlling sensor are approximately 
1.6 and 3.6 prad rms for the X and Y axes, respectively. 

Noise from the Atlas missile gyros is characterized by numerous 
discrete frequency components on top of a generally white noise background. 
The discrete power spikes are peculiar to ball bearing gyros. The dis- 
crete components contribute over 85 percent of the overall gyro rms 
noise level of 6.2 prad in the X axis, and 12.4 prad in the Y 
axis. Several of the discrete components arise from sources related to 
the wheel spin frequency. In particular, a strong 71 Hz component 
appears to be related to the rotation frequency of the rotor spin bearing 
retainers. 

The primary autocollimator noise spectrum is essentially white, 
and is due to shot noise created in the photodetector. The rms noise 
equivalent angle in the frequency band 0.1 to 100 Hz is 0.5 prad. 

In contrast to the white primary autocollimator noise, the azimuth ref- 
erence autocollimator has primarily l/f noise. The rms noise equivalent 
angle is 2.8 prad. 

B* Application to Relativity Experiment 

It is likely that a test of the actual relativity experiment sensors 
will require active pointing of the experiment package to keep signals 
within the range of intended operation. Consequently, spectra obtained 
will be dependent on pointing loop dynamics, and proper interpretation 
will require a system model similar to the one developed for the Fixed 
Base Simulator. The same computer programs developed to compute the 
spectral functions of simulator signals will be directly applicable to 
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the new system. The present data acquisition system is cumbersome, and 
it is recommended that a direct digital acquisition system be considered. 

With the increased precision of the relativity sensors, it is likely 
that base motion will contribute a much larger portion of the sensor 
output signals. The relative amount of base motion can be determined 
by using the coherence function 

The large spikes observed in the simulator gyro spectra, caused by 
bearing irregularities, are not likely to be present in the relativity 
gyro signals. However, sinusoidal noise or systematic components produced 
either in the pickoff or the rotor will be detected. Care should be taken 
in determining the amplitude of such components from spectral estimates 
due to the effects of data windows, prewhitening, and other computational 
considerations discussed in Appendix A. 
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Appendix A 

COMPUTING SPECTRAL FUNCTIONS 


Straightforward application of the formulae of Chapter IV to esti- 
mate discrete spectral functions is generally neither efficient nor reli- 
able in producing useful results. Several considerations must be made 

c 

to minimize the computation required and to account for the finite length, 
relative power distribution, and randomness of the data. These consider- 
ations are outlined below. Three references [B5], [B2], and [K2 ] in par- 
ticular are recommended for further reading in interpreting spectral 
plots and in avoiding numerous computational pitfalls. 

1. Fast Fourier Transform 

An initial inspection of the discrete Fourier transform formula (4.5) 

indicates that the number of computer operations required is proportional 
2 

to N . However, as discussed in [C3] several algorithms have been devel- 
oped which require a number proportional to N log N. The ratio of the 
number of operations or equivalently of the time required is log N/N. 

Dramatic reductions in computation time can therefore be achieved when 
N is large. As an example, when N = 1024, the new algorithms require 
only about 1/100 the time. 

The most widely known and used algorithm is the Fast Fourier Trans- 
form (FFT) developed by Cooley and Tukey [C4], A considerable body of 
literature exists on the FFT, including the entire June 1967 special issue 
of the IEEE Transactions of Audio and Electroacoustics . As a consequence, 
a description of the algorithm will not be given here. Intuitive devel- 
opments and explanations of the FFT are given in [ B8] , [Gl], [D3], and 
[ G3] . 

2. Leakage 

A consequence of finiteness of data is the loss of frequency resolu- 
tion in the spectral functions. Consider a continuous finite length sig- 
nal x ( t ) . This signal can be viewed as an infinitely long process mul- 
tiplied by a time function, called a data window, which has unit value 
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over the finite signal length and is zero elsewhere 

x f (t) = d(t) x(t) 

where 

x(t) = infinitely long signal 

/0 t < 0 

d(t) = 1 0 < t < T 

*0 t > T 

x f (t) = finite length record 

T = record length 

The effect of this time multiplication is the equivalent in the frequency 
domain of convolving the Fourier transform of x(t) with the Fourier 
transform of d(t). The latter transform is called a spectral window. 

The spectral window association with the rectangular data window d(t) 
is easily shown to be of the form sin (af)/(af) where a is a constant. 
Because of the convolution, spectral estimates at a point are spread over 
neighboring frequencies. Using the rectangular data window, the contri- 
bution of a signal component at frequency f = f^ + Af to the power spec- 
tral estimate at f n falls off as l/|Af| 2 . Conversely, a signal at f^ 
affects neighboring estimates as l/|Af| 2 . This phenomenon is known as 
leakage and can seriously affect the validity of power spectrum estimates. 

Leakage also occurs in transforming discrete signals. Here, however, 
a signal component exactly at one of the discrete frequencies of the Dis- 
crete Fourier Transform (DFT) does not affect the estimate at any other 
frequency. On the other hand, signal components not at a discrete fre- 
quency of the DFT display leakage just as in the continuous case. Since 
most of the signals of interest do not fall at discrete frequencies of 
the DFT, it is necessary to reduce the effects of leakage by applying a 

different data window which will cause leakage to fall off at a rate 

2 

greater than l/|Af| without broadening the main lobe of the spectral 
window excessively. 

Blackman and Tukey, in their comprehensive book on computation of 
power spectra [B5], analyze several different windows. Selections widely 
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used are the "Hanning" window 


d (t) 
m 


/ 0.5(1 - cos 2flt/T) 
VO 


0 < t < T 
otherwise 


and the " Hamming" 


window 


d „(t) 


0,54 - 0.46 cos 2itt/T 0 < t < T 
0 otherwise 


Both data windows yield a spectral window with a narrow main lobe and 

small (1$> - 2 <f) side lobes. However, harming results in leakage falling 

off as l/|Af[ , for TAf » 1, whereas hamming results in leakage atten- 

2 

uation of only approximately l/[Af| for TAf » 1. The advantage of 
hamming is that the highest spectral window side lobe is one third as high 
as that of the hanning window. For spectra which do not rapidly change 
with frequency, that is where leakage is not of major concern, the lower 
peak side lobes of the hamming window can be more desirable than the in- 
creased leakage rolloff obtained from hanning. 

A third window, called the Bartlett window, is described by 

d b (t) = 1 - 1 0 < t < T 

= 0 otherwise 


The main advantage of this window is the ease of computation, since it 

4 

requires no trigonometric functions. Leakage falls off as 1/ | Af | . 

Choice of an appropriate data window depends upon the form of the 
power spectrum. It should be noted, however, that any of the above win- 
dows will work satisfactorily in most applications where the spectrum 
changes slowly with frequency. 

Multiplying the data by a window decreases the power in the esti- 
mate. Therefore, to maintain the proper amplitude in the power spectrum 


165 



estimate, it is necessary to normalize the estimate by the mean of the 
squared window. That is 


where 


yN Af 

P (n) = — 
XX u 


N-l 

I 

k=sO 


x(k) d(k) W nk 


N-l 

j=0 


d 2 (j) 


3. Prewhitening 

As a result of leakage, significant errors can be introduced in the 
estimate when power is concentrated at a particular frequency. For ex- 
ample, test data almost always has some very low frequency noise including 
DC. This noise may come from instrument and data acquisition equipment 
biases and drifts or from analog to digital converter offsets. Due to 
leakage, the true values of the power spectrum at the low frequencies may 
be completely overcome by the relatively large side lobes of the spectrum 
at or near DC. it is desirable, therefore, to remove large amounts of 
noise at particular frequencies before computing the power spectrum. This 
procedure is known as prewhitening, since the resulting spectrum is more 
nearly white. 

The prewhitening technique used in analyzing the Fixed Base Simulator 
consists of fitting a straight line and/or a sine wave to the data using 
unweighted least squares. The resulting fits are then subtracted from the 
data, and the sample mean and variance of the residual data computed. Since 
bad data points sometimes occur, all points greater than 4 a, where a is 
the sample standard deviation, are output and then set equal to the mean of 
the residual data. This effectively eliminates bad data without signifi- 
cantly affecting the power spectrum. if any points are found to exceed 
4a, the entire prewhitening procedure is repeated with the bad points elim- 
inated. 
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4. 


Randomness 


The effects of data randomness must also be considered in addition 
to the frequency distortion due to finite data and strong frequency com- 
ponents. In Chapter IV, it was pointed out that forming spectral esti- 
mates from increasing amounts of data does not necessarily improve the 
estimate. In other words, the estimate of the spectrum may not converge 
to the true value no matter how much data is processed. 

Take, for example, the case where the spectrum of the signal is 
white; i.e., it is constant for all frequencies. The estimate of the 
power at zero frequency is from (4.17), 


P(0) 


At 

N 


N— 1 

I 


k=0 


x(k) 


2 


The x(k) are assumed independent and Gaussian with zero mean and vari- 
2 

ance a . Let At = 1, for convenience. The expected value of the esti- 
mate is 


Ef(0» - ^ E 


N-l 

I 

k=0 


2 


x(k) 


2 

a 


The variance of P(0) is 


or 


Var P(0) = E{P(0) - E[P(0)] 


} 


2 


Var P(0) 



N-l 


I 

k=0 


x(k) 


4 


4 

a . 


Application of the moment generating function of a normal process yields 

Var P(0) = 2 q 4 
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The ratio of the standard deviation of P(0) to its expected value is, 
then, 


x/Var P(0) _ i— 

E[P(0)] ' 

in other words, the rms error in estimating the actual P(0) is greater 
than the expected value by a factor of sfe, independent of record length 
N. Obviously, means of reducing such large uncertainties in the power 
spectrum have to be employed. 

Although many approaches to the problem of reducing statistical 
variability have been discussed [B5], [B4], and [R2], the one selec- 
ted for use is that published by Welch [W2], This technique divides the 
entire data record into shorter, possibly overlapping subrecords, applies 
a data window to reduce leakage, computes a power spectrum for each sub- 
record, and averages the results point by point. Another approach is to 
convolve the spectrum obtained from one record with a spectral window. 

The convolution effectively forms a weighted average of the spectral esti- 
mates near the frequency of interest. The smoothed results obtained are 
achieved in this manner by a corresponding loss in frequency resolution. 

The approach of Welch was selected since it most efficiently uses 
all of the data and yields finer resolution than the frequency convolu- 
tion method. Improved resolution results not because of properties of 
the data, but because better use could be made of the limited storage 
capability of the XDS Sigma-5 computer on which the computations were 
made. The maximum number of positive frequency values which can be com- 
puted at one time on the Sigma-5 is 2049. Therefore, if a sufficient 
number of spectra of subrecords are averaged to achieve statistical sta- 
bility, the resolution is f^/2048, where f^ is the sampling frequency. 

To achieve the same resolution using the frequency convolution method, 
a spectrum containing M X 2048 points would have to be initially computed, 
where M is the number of subrecord spectra computed. This is because 
the resolution of the average-over-frequency method is reduced by a fac- 
tor about equal to the number of spectral estimates averaged. Since the 
maximum length of spectrum is specified by computer constraints, the 
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resolution available using average-over-frequency methods is MAf, where 
Af is the resolution of the average- over- time method. 

How averaging several spectral estimates improves stability can be 
analyzed by noting that the spectral estimates P(n) are normally dis- 
tributed for normally distributed noise data. The average of several 
P(n) therefore has a chi-square probability distribution, where 

2E 2 P(n) 

Var P(n) 


is the equivalent number of degrees of freedom. Averaging over M non- 
overlapping subrecords gives k = 2M, or expressed differently, 


Var P(n) 


E 2 P(n) 
M 


For example, averaging 25 independent spectral estimates will yield a 
value within 20 $ of the true value with 67$, confidence. This relation- 
ship holds only for normally distributed noise, but can be used as a 
rule of thumb for other distributions. Using this expression a tradeoff 
between relative stability and computation time can be made. 

5. Sampl ing 

Since the pow.er spectra are computed digitally, it is necessary to 
digitize analog test data. The FFT algorithm assumes that the data points 
are equally spaced in time, and so it is convenient to sample at a fixed 
frequency. The act of sampling, however, can make spectral estimates 
meaningless unless noise at frequencies greater than half the sampling 
frequency is first rejected. 

It is well known that the power spectrum of an ideally sampled sig- 
nal is 


P*(f) = P(f - kf ) 
s 

where P(f) is the power spectrum of the original continuous signal, 

and f is the sampling frequency. This implies that high frequency 
s 
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noise in P(f) can appear as low freqeuncy noise in the sampled spec- 
trum P* (f) . If, for example, the signal is white noise bandlimited 
at ten times the sampling frequency, the resulting sampled spectrum at 
low frequencies is ten times greater than the actual. 

Distortion of this nature, which is known as aliasing, is not present 

when P(f) =0, | f I > f /2. For then P* (f) = P(f) , -f /2 < f < f /2. 

s s s 

Passing the signal through a low-pass analog filter with bandwidth less 

than f /2 therefore effectively eliminates aliasing, 
s 

The three figures A. 1, A. 2 and A. 3, demonstrate the effectiveness 
of filtering. A 10 Hz square wave from a low-frequency function gen- 
erator was passed directly into an analog-to-digi tal converter (ADC) 
sampling at 409.6 Hz. Ideally, the power spectrum should be composed of 

all the odd harmonics of the fundamental frequency with amplitudes falling 
2 

off as 1/f . This can be seen in Fig. A. 1 which is the spectrum of a com- 
puter generated 5. 12 Hz square wave. In practice even harmonics and noise 
will appear due to non-ideal aspects of the ADC, function generator, and 
signal transmission. The power spectrum of the unfiltered sampled signal 
is shown in Fig. A. 2. The original analog signal was then passed through 
a 100 Hz low-pass filter before analog-to-digital conversion. The result- 
ing spectrum is shown in Fig. A. 3. Comparison of Figs. A. 2 and A. 3 shows 
the extensiveness of aliasing present in Fig. A. 2. Except for the small 
even harmonics, the spectrum in A. 3 is nearly ideal up to 100 Hz. 

0. Digital Filtering 

Once a signal has been sampled, there is no way to remove aliasing 
except by resampling the analog signal either at a higher frequency or 
after low-pass filtering. However, a low-passed analog signal sampled 
at one frequency can be sampled at a lower rate by first eliminating higher - 
frequencies in the digitized information using a low-pass digital filter. 
The sampling rate is reduced by skipping points or averaging several points 
of the digitally filtered signal. Aliasing can still be introduced if the 
new sampling rate is less than twice the highest frequency present in the 
digitally filtered signal. Digital filtering is especially useful in ob- 
taining low-frequency spectral information when the original analog signal 
is sampled at a high frequency. 
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Fig. A. 1. POWER SPECTRUM OF A COMPUTER GENERATED 5.12 Hz SQUARE WAVE 
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7. 


Quantization 


Another important aspect of sampling is the necessity of quantizing 

the data so that it can be represented as a digital word. in the A/D 

conversion process, the voltage level of the signal must be represented 
N 

by one of 2 possible numbers, where N is the number of bits per word. 

A roundoff error is therefore introduced into the digital data. Widrow 

[W3] has shown that this quantization error can be modeled as additive 

white noise, uniformly distributed between -q/2 and q/2 with zero 

2 

mean and variance q /12, where q is the quantization level. Using 
this model, the contribution of quantization error to the power spectrum 
is 


P(n) = q 2 At/12 

Quantization error can be minimized by scaling the maximum signal level 
to the largest value which can be digitized. 

8. Graphical Representation 

A log-log plot of the spectral functions has been found to be the 
most useful presentation of the computed results. Most signal analysts 
are familiar with this form and can readily interpret the slopes of 
straight line segments and break frequencies in terms of their underlying 
physical mechanisms. The logarithmic scales also allow a wider range of 
data to be presented on the same plot. 

One effect of logarithmic frequency scaling is the graphical com- 
pression of the equally spaced spectral estimates with increasing fre- 
quency. Since equal frequency resolution is normally not required across 
the entire spectrum, it has been suggested [K2] that only the average of 
several points in a frequency band be displayed, with more points averaged 
with increasing frequency. One convenient way to maintain an approximately 
constant fractional frequency resolution is to break up the spectrum into 
bands each twice as large as the next lower band and display only the av- 
erage of the points in each band. Since more points are averaged at higher 
frequencies, greater smoothness results. A negative aspect is that a spike 
in the spectrum due to a discrete component can be lost due to the averaging. 
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As a compromise, the data is averaged as above, but any points greater 
than three times the average value are also plotted. 
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Appendix B 


SCHEMATIC DIAGRAMS OF THE FIXED BASE SIMULATOR 

Schematics of the electronic circuits developed for the Fixed 
Base Simulator are included in the following pages. Indicated card 
numbers correspond to those described in Chapter II and to the gyro and 
autocollimator preamplifiers mounted on the platform. Diagrams of the 
mode-switching relays and the meter panel are also included. 
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CARDS 1 & 3: GYRO TORQUER AMPLIFIER 
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CARD 2: TORQUER AMPLIFIER DECOUPLING NETWORK 
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CARD 4: ±6.2 VDC REGULATOR AND SENSOR BIASING NETWORK 
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CARD 4: GYRO TORQUE REBALANCE (CAGING) AMPLIFIERS 
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CARD 5: MODE CONTROL LOGIC AND DRIVERS 
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CARD 7: GYRO HEATER CONTROL 
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CARD 8: AUTOCOLLIMATOR ELECTRONICS 
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CARD 10: A/C INSTRUMENTATION FILTERS; ±12 VDC POWER SUPPLY 
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CARD 21 



GYRO PREAMPLIFIERS 
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CARD 22: GYRO PREAMPLIFIER FILTER 
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METER PANEL 
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MODE SWITCHING RELAYS - RELAY PANEL 
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